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OBJECTIVES:  To  characterize  microbial  communities  in  biofilms 
lining  marine  infaunal  burrows;  to  isolate  new  reductively 
dechlorinating  marine  bacteria  from  burrow  lining  biofilms; 
to  determine  microscale  placement  of  reductive  dechlorinators 
in  biofilms. 

APPROACH:  Enrichment  cultures  providing  a  range  of  growth 
conditions  were  inoculated  with  burrow  and  tube  lining 
biofilm  material  collected  in  the  field  from  bromoaromatic 
producing  Balanoglossus  aurantiacus)  and  non-producing 
{Kinbergonuphis  jenneri)  infauna.  These  enrichments  were 
screened  I  for  dechlorination  of  2,4,6-trichlorophenol  and  were 
sources  for  isolation  of  dechlorinating  pure  cultures.  The 
placement  of  microorganisms  in  tube  lining  biofilms  has  also 
been  examined.  Since  various  infaunal  burrows  and  tubes 
support  remarkably  similar  microbial  communities,  we  have 
focused  on  a  single  model  tube  system  constructed  by  the 
onuphid  polychaete  Diopatra  cuprea.  We  used  confocal  scanning 
laser  microscopy  to  visualize  the  microbiota  in  Diopatra  tube 
biofilms  and  to  determine  which  of  these  organisms  are 
potentially  active.  A  neural  computing  approach  to  analysis 
of  phospholipid  fatty  acid  profiles  has  been  developed.  Pure 
cultures  of  reductively  debrominating  marine  bacteria  have 
been  characterized  and  the  reductive  debrominase  they 
expressed  examined. 

ACCOMPLISHMENTS :  Several  enrichment  culture  series,  prepared 
in  1997  and  1998  from  field-collected  samples  of  K.  jenneri 
tubes  and  B.  aurantiacus  burrow  linings,  completely 
mineralized  2,4,6-trichlorophenol  and  yielded  several 
anaerobic  pure  culture  isolates  that  reductively 
dechlorinated  this  compound.  These  isolates  lost  the 
capacity  to  dechlorinate  after  several  transfers.  We  tried  a 
variety  of  vitamin  supplements  and  basal  medium  recipes,  but 
were  unable  to  restore  dechlorination  activity  to  these 
isolates,  New  enrichments  have  been  prepared  using  more 
restrictive  medium  recipes.  These  enrichments  have  performed 
well  through  initial  transfers  and  have  yielded  one 
dechlorinating  pure'  culture  to  date.  This  organism  is 
currently  being  characterized  and  our  isolation  efforts  are 
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continuing.  We  have  completed  the  characterization  of  two 
reductively  debrominating  strains  of  Propionigenium  marls. 
These  cultures  convert  mM  levels  of  2,4,6-tribromophenol  to 
2-  or  4-bromophenol  and  were  isolated  from  the  burrows  of  the 
hemichordate  B.  aurantiacus  and  the  capitellid  polychaete 
Notomastus  lobatus.  These  strains  differ  in  a  number  of 
structural  and  physiological  characteristics  from  the  P. 
marls  type  strain,  particularly  in  the  rapid  debromination  of 
2,4,6-tribromophenol  by  our  strains  and  the  near  absence  of 
debromination  by  the  type  strain.  Although  these  strains  are 
distinguishable  on  the  basis  of  several  key  characteristics, 
the  16S  rRNA  sequence  homology  among  them  is  greater  than 
99%.  This  may  indicate  rapid  evolution  of  burrow-associated 
P.  marls  strains.  The  reductive  debrominase  produced  by 
these  strains  is  cytoplasmic  and  can  be  reduced  by  both  NADH 
and  NADPH.  There  is  no  evidence  for  P.  marls  respiration 
using  bromophenols  as  terminal  electron  acceptors. 

We  had  previously  examined  the  biofilms  lining  burrows  and 
tubes  of  several  polychaetes  and  hemichordates  using 
phospholipid  fatty  acid  (PLFA)  analysis  and  have  now  employed 
a  neural  computing  approach  to  resolve  differences  among  PLFA 
profiles.  This  analytical  method  confirmed  the  similarity  of 
bacterial  communities  in  tubes  and  burrows  of  a  variety  of 
animals,  allowing  us  to  identify  and  focus  on  a  single  model 
infaunal  tube  system.  We  have  chosen  tubes  constructed  by  a 
ubiquitous  onuphid  polychaete,  Dlopatra  cuprea,  as  a  model 
experimental  system.  This  animal  is  found  across  a  range  of 
sediment  types  and  can  be  collected  from  chronically 
contaminated  sites.  Tubes  built  by  a  related  onuphid, 
Klnbergonuphls  jennerl,  have  also  been  used  in  some 
experiments. 

We  have  developed  mounting  and  fluorescent  staining  methods 
for  confocal  scanning  laser  microscope  examination  of  the 
biofilms  lining  Dlopatra  cuprea  tubes.  The  average  thickness 
of  these  biofilms  was  about  60  |J.m.  No  channels  or  channel¬ 
like  structures,  which  are  common  in  most  laboratory- 
cultivated  biofilms,  were  observed.  Bacteria  occurred  as 
both  single  cells  and  microcolonies  throughout  the  biofilm 
material.  Single  cell  numbers  averaged  about  5.6  x  10®  ‘cm"® 
and  the  majority  of  these  cells  (ave.  68%)  were  intact,  had 
membrane  potential  (as  found  using  Syto  9/propidium  iodide 
staining) ,  and  presumably  capable  of  metabolic  activity.  No 
trends  in  active  cell  numbers  with  increasing  depth  in  the 
biofilm  were  found,  indicating  potential  for  activity 
throughout  the  structure.  Microcolony  sizes  were  highly 
variable,  ranging  from  16  |J.m®  to  2.5  x  10®  (Xm®.  Single  cell 
size  (ave.  =  0.085  jxm®)  decreased  and  became  more  uniform 
with  depth  in  the  biofilm.  In  situ  hybridization  strategies 
for  localization  of  reductively  dechlorinating  and  other 
anaerobic  bacteria  have  been  developed. 


CONCLUSIONS ;  Biofilms  lining  burrows  and  tubes  constructed  by 
marine  infauna  contain  a  diverse  and  abundant  microbiota, 
occurring  as  single  cells  and  in  microcolonial  formations . 

The  majority  of  these  organisms  are  at  least  potentially 
active.  The  burrow  biofilm  microbiota  includes  reductively 
dehalogenating  bacteria,  which  can  be  readily  enriched  from 
burrow  structures  of  several  different  infaunal  species.  At 
least  some  of  these  bacteria  are  amenable  to  isolation  and 
can  be  cultivated  in  the  laboratory.  We  have  enriched 
reductively  dechlorinating  bacteria  from  burrows  produced  by 
both  halometabolite  producing  and  non-producing  marine 
infauna.  It  is  clear  from  these  results  that  marine  - 

dechlorinating  bacteria  are  common  in  infaunal  burrow 
structures  relative  to  bulk  sediments.  It  is  also  clear  from 
our  microscopic  examination  of  tube  lining  biofilms  that  the 
potential  for  interspecific  interactions  among  dehalogenating 
bacteria  and  other  types  of  organisms  exists.  These 
interactions  may  include  consortial  metabolism  of  chlorinated 
aromatic  compounds  by  combinations  of  reductive 
dechlorinators  and  aromatic  ring  cleaving  organisms. 

SIGNIFICANCE;  This  study  has  revealed  and  characterized  an 
important  and  previously  unrecognized  niche  for  reductively 
dehalogenating  bacteria  in  marine  sediments.  Infaunal  burrow 
structures  provide  a  source  for  new  organisms  with  potential 
application  in  bioremediation  processes. 
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ABSTRACT 


The  microbial  community  composition  of  surface  and  subsurface  marine  sediments,  and 
sediments  lining  burrows  of  marine  worms  from  the  North  Inlet  estuary  were  analyzed  by 
comparing  ester-linked  phospholipid  fatty  acid  (PLFA)  profiles  using  neural  network  (NN)  and 
cluster  analysis.  PLFA  profiles  were  converted  to  binary  format  for  transmission  to  a  back 
propagating  NN.  The  NNs  were  trained  to  relate  binary  input  to  sample  type  (e.g.,  surface, 
subsurface,  and  burrow  lining)  and  worm  species  (e.g.,  Notomastus  lobatus,  Balanoglossus 
aurantiacus  and  Branchyoasychusamericana).  The  optimal  architecture  for  the  NN  was  found  to 
be  540-input,  16-hidden  and  4-output  neurons.  The  optimal  amount  of  training  was  found  to  be 
approximately  150  epochs  (or  to  an  error  tolerance  of  0.01).  Results  from  cross-validation 
revealed  that  NNs  were  able  to  correctly  discriminate  microbial  community  profiles  primarily  by 
sample  typie  (86.4%)  rather  than  worm  spiecies  (66. 1%).  Dendrograms  produced  by  using 
Euclidean  distance  and  average  linkage  clustering  methods  on  data  from  five  independently  trained 
NNs  yielded  the  following  results,  (i)  PLFA  profiles  from  microbial  communities  in  the  surface 
sediments  were  characterized  by  high  levels  of  long-chain  (>  20  carbon  units)  polyunsaturated 
PLFAs  and  low  levels  of  methyl-branched  saturated  PLFAs.  (ii)  Communities  associated  with 
subsurface  sediments  were  characterized  by  low  levels  of  polyunsaturated  PLFAs  and  high  levels 
of  methyl-branched  PLFAs.  (iii)  Sediments  lining  burrows  of  marine  worms  had  PLFA  profiles 
that  were  intermediate  between  these  two  extremes  presumably  reflecting  the  oxic/anoxic  interface. 
The  utility  of  NNs  for  interpreting  complex  PLFA  profiles  of  microbial  communities  was 
demonstrated. 
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INTRODUCTION 

Natural  microbial  communities  are  complex  assemblages  of  organisms  composed  of  a 
variety  of  different  physiological  groups  of  bacteria,  archaea,  tind  microeucaryotes,  including 
fungi,  cilgae  and  protozoa.  Recent  molecular  biological  studies  have  shown  that  the  assemblage  of 
organisms  present  in  an  environmental  sample  can  contain  thousands  of  distinct  biotypes  of 
bacteria  (45, 46)  and  that  many  of  these  organisms  represent  undescribed  microbial  lineages,  some 
of  which  lack  even  a  single  representative  species  available  in  pure  culture  (3).  This  complexity  is 
further  compounded  by  the  high  spatial  and  temporal  variability  of  species  represented  in  many 
natural  communities  (1,6),  and  by  the  still  open  question  of  what  proportion  of  the 
microscopically  observable  assemblage  of  organisms  is  viable  or  even  metabolically  functional  (8, 
29, 53).  Characterization  of  natural  microbial  communities  is  clearly  a  difficult  problem.  Many 
molecular  biological  methods  for  community  analysis  are  either  very  time  consuming  such  as 
construction  and  screening  of  clonal  libraries  (3, 11, 20),  or  yield  an  index  of  community 
composition  that  does  not  reveal  organism  taxonomic  affiliations  without  additional  analysis,  such 
as  denaturing  gradient  gel  electrophoresis  (16, 34, 41).  There  is  a  clear  need  for  methods  than  can 
be  used  to  rapidly  and  accurately  examine  microbial  community  composition.  Ideally,  such 
methods  would  include  all  (or  at  least  most)  microbial  domains  and  produce  a  standardizeable 
community  profile,  allowing  comparisons  across  different  sample  types  and  facilitating 
development  of  broad  hypotheses  concerning  microbial  community  dynamics. 

One  approach  to  microbial  community  characterization  that  has  been  successful  in  a  wide 
variety  of  applications  (e.g.,  4,  5,  9,  18,  21, 35, 37,  43),  is  phospholipid  fatty  acid  (PLFA) 
analysis  (17).  Phospholipids  are  structural  components  of  all  biological  membranes.  These 
compounds  have  no  storage  function  and  thus  represent  a  consistent  fraction  of  cell  mass.  They 
also  degrade  quickly  upon  an  organism's  death  and  current  extraction  and  derivatization  methods 
permit  recovery  of  PLFAs  exclusively  from  living  organisms.  Extraction  and  subsequent  analysis 
by  gas  chromatography/mass  spectroscopy  provides  precise  resolution,  sensitive  detection,  and 
accurate  quantification  of  a  broad  array  of  PLFAs.  Analyses  can  be  performed  rapidly  and 
effficiently  and  each  yields  a  profile  composed  of  numerous  PLFAs  defined  on  the  basis  of 
compound  structure  and  the  quantity  of  each  compound  present  in  the  sample.  Different  PLFA 
average  chainlengths  and  degrees  of  unsaturation  are  produced  by  bacteria  and  eucarya,  allowing 
both  domains  to  be  examined  through  the  same  analysis.  The  PLFA  profile  obtained  thus 
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constitutes  a  'fingerprint'  of  the  living  bacterial  and  microeucaryotic  microbial  community  and 
reflects  its  species  composition.  Comparisons  across  sample  types  are  simplified  by  normalizing 
the  quantity  of  each  PLFA  to  the  total  recovered,  expressing  the  result  in  terms  of  mol%  PLFA. 
One  drawback  of  the  PLFA  analysis  is  that  archaea  are  not  represented  in  the  analysis  since  their 
membrane  lipids  employ  ether  rather  than  ester  bonds.  Nonetheless,  PLFA  profiles  can  be  used  to 
characterize  microbial  community  compositions  and  provide  useful  information  on  the  dynamics  of 
viable  bacteria  and  microeucaryotes. 

Comparing  PLFA  profiles  of  microbial  communities  is  often  complicated  by  variability  in 
the  types  and  quantities  of  compounds  present  in  replicate  samples  from  a  given  sample  site.  This, 
together  with  the  complexity  of  the  PLFA  profiles,  constitutes  substantial  challenges  to  our  ability 
to  quantitatively  compare  microbial  community  compositions  using  this  approach.  Previous 
studies  have  employed  cluster  analysis  and  principal  components  analysis  (PCA)  to  group  like 
samples  (43)  and  to  determine  which  PLFAs  are  most  diagnostic  for  the  groups  defined  (22, 25, 
42).  However,  these  methods  do  not  permit  details  of  the  profiles  to  be  conveniently  examined.  A 
method  that  provides  the  advantages  of  cluster  analysis  and  PCA  while  facilitating  more  detailed 
profile  examination  would  be  very  advantageous. 

Artificial  neural  networks  (NNs)  provide  an  useful  tool  for  recognizing  patterns  in  complex 
data  sets  such  as  those  associated  with  PLFA  'fingerprints'.  NNs  are  constructed  by  using 
computer  software  and  consist  of  layers  of  neurons  that  make  independent  computations  and  pass 
on  their  outputs  to  other  neurons  (36).  Each  neuron  in  a  layer  is  connected  to  neurons  in  the  next 
layer  so  that  output  of  each  neuron  affects  the  activation  of  all  neurons  to  which  it  is  connected. 
Neurons  are  adaptable,  and  through  the  process  of  learning  from  examples,  store  knowledge  and 
make  it  available  for  use  (2).  In  a  training  technique  called  error  back-propagation  (38),  a  pattern  is 
presented  to  an  input  layer  of  a  network,  and  the  network  produces  output  based  on  the  sum  of  the 
weighted  inputs.  When  the  pattern  of  the  output  layer  is  compared  to  target  values,  the  errors 
between  them  are  computed.  An  error  function  is  used  to  readjust  the  weights  of  each  neuron. 

This  iterative  process  continues  until  the  errors  converge  to  a  value  below  a  designated  tolerance 
level  (38).  The  adjusted  weights  of  a  trained  network  can  be  used  to  recognize  patterns  such  as 
those  in  PLFA  profiles  of  microbial  communities  and  to  provide  functional  relations  between 
profiles  and  target  values. 
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There  is  a  paucity  of  information  concerning  the  effects  of  NN  architecture  and  training  on 
NN  convergence  and  prediction.  For  the  non-specialist,  selecting  the  NN  architecture  and 
determining  the  amount  of  training  needed  for  optimal  performance  of  NNs  has  not  been  clearly 
established.  For  example,  several  studies  have  pointed  out  that  having  too  many  neurons  in  the 
hidden  layer  may  cause  unnecessary  memorization  of  the  training  set  (27, 28, 39, 51).  As  a 
consequence,  the  NN  loses  the  capacity  to  generalize.  Having  too  few  neurons  in  the  hidden  layer 
results  in  non-convergence  of  the  training  set  to  the  desired  error  tolerance.  Overtraining  NNs  to 
low  error  tolerance  values  has  similar  consequences. 

The  focus  of  this  study  was  to  compare  PLFA  profiles  of  microbial  communities  in  surface 
and  subsurface  nearshore  sediments  and  in  the  linings  of  burrows  of  three  species  of  marine 
worms  by  using  a  back-propagating  neural  network,  and  to  outline  the  procedures  needed  to 
optimize  NN  performance.  Optimization  was  achieved  by  varying  the  architecture  and  the  extent  of 
training  and  determining  the  effects  of  these  factors  on  NN  prediction.  Here,  we  report  on  two 
different  neural  computing  approaches  for  comparing  PLFA  profiles:  cross-validation  and  cluster 
analysis  using  data  from  trained  NN. 

MATERIALS  AND  METHODS 

Source  of  PLFA  profiles  analyzed.  PLFA  data  used  in  our  analyses  were  from  a  previous 
study  of  marine  sediment  microbiota  by  Steward  et  al.  (42).  This  study  should  be  consulted  for 
sample  site  descriptions  and  significant  features  of  the  internal  burrows  sampled.  Briefly, 
sediment  samples  collected  for  that  study  were  ten  matched  sets  for  each  of  three  worm  species, 
with  each  set  including:  1  worm  burrow  sediment  sample,  1  surface  sediment  sample,  and  1 
subsurface  sediment  sample.  Surface  sediments  were  collected  from  approximately  the  top 
centimeter  of  sediment  near  but  not  adjacent  to  the  worm  burrows.  Subsurface  samples  were  taken 
approximately  10-20  cm  deep,  near  but  not  adjacent  to  worm  burrows.  Sediments  were  scraped 
from  the  linings  of  burrows  built  by  the  capitellid  polychaete  Notomastus  lobatus  and  the 
IxemxchorddXe  Balanoglossiisaurantiacus  using  sterile  spatulas.  Segments  of  tubes  constructed  by 
the  maldanid  polychaete  Branchyoasychus americam  were  recovered  intact.  Large  macrofauna 
were  excluded  from  all  samples.  The  samples  were  placed  in  sterile  Whirlpack  bags  and 
immediately  frozen  on  dry  ice,  then  transferred  to  -80®C  pending  PLFA  extraction  and  analysis  as 
described  by  Steward  et  al.  (42).  This  sample  set  included  90  total  samples,  each  yielding  a  profile 
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of  about  60  different  fatty  acids.  Forty-two  PLFAs  were  present  at  levels  >  0.01  niol%  of  most 
profiles  accounted  for  93-97%  of  total  PLFA  in  each  sample  (42).  The  balance  represented  low 
abundance  compounds. 

Data  for  the  NN.  PLFA  data  used  in  this  analysis  also  included  control  data  which  consisted  of 
23  randomly  generated  PLFA  profiles.  The  limits  of  the  random  numbers  were  set  by  the 
maximum  and  minimum  values  of  the  PLFA  data. 

Input  data  for  training  the  NN  were  generated  by  converting  each  of  the  1 12  PLFA 
components  to  9-digit  binary  numbers.  This  conversion  set  the  maximum  and  minimium  PLFA 
values  to  between  '(XXXXXXXX)'  and  '111111111',  respectively.  For  example,  if  the  quantity  of  the 
PLFA  16:0  was  2.71,  the  converted  binary  designation  would  be  '0001 1(XX)0'.  Input  data  for  an 
entire  microbial  community  in  a  sample  consisted  of  60  9-digit  binary  numbers  ordered  by  their 
detection  by  the  mass  spectrophotometer. 

Output  data  for  training  the  NN  related  input  data  to  sample  type  and  sample  group  (by 
worm  species)  and  was  coded  as  a  4-digit  binary  number.  The  first  2-digits  refer  to  the  worm 
species  and  the  next  2-digits  refer  to  sample  type.  The  worm  species  were  coded  '00',  '10',  '01' 
and  'll'  corresponding  to  the  marine  worms B. aurantiacus,  B. americana,  and N.  lobatus and 
to  randomly-generated  data,  respectively.  The  sample  types  were  coded  '01',  '00',  '10'  and  'll' 
and  corresponding  to  surface,  burrow,  and  subsurface  samples,  and  randomly-generated  data, 
respectively. 

Input  and  output  data  used  for  training  the  NN  consisted  of  the  60  9-digit  binary  numbers 
and  the  4-digit  binary  output  data.  Data  used  for  'testing'  the  NN  consisted  of  the  60  9-digit  binary 
numbers  only.  To  conform  with  the  formating  requirements  of  the  NN  program,  spaces  were 
inserted  between  each  digit. 

Cross-validation  scheme.  To  evaluate  the  predictive  power  of  the  neural  network,  we 
employed  the  following  cross-validation  scheme:  the  order  of  the  data  were  randomized,  75%  of 
the  data  were  used  to  train  to  an  average  error  per  sample  of  0.01,  and  the  remaining  25%  were 
used  to  validate  the  NN  performance.  Correct  and  incorrect  classifications  were  recorded.  The 
above  scheme  was  repeated  10  times  for  each  cross-validation. 

Optimization  of  the  NN.  The  effects  of  architecture  on  the  predictive  accuracy  of  the  NN  were 
determined  by  independently  training  the  NNs  using  4, 8, 16, 64,  or  270  neurons  in  the  hidden 
layer.  Unless  otherwise  specified,  all  NNs  were  trained  to  an  eitor  tolerance  of  0.01  using  a 
learning  rate  of  0. 1 ,  a  momentum  parameter  of  0,  and  a  noise  factor  of  0.  Duplicate  cross- 
validation  experiments  were  used  to  evaluate  the  ability  of  trained  NNs  to  relate  PLFA  profiles  to 
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sample  types  and  worm  species. 

To  determine  the  effects  of  error  tolerance  level  on  the  predictive  accuracy  of  the  NN,  error 
tolerance  values  of  0.5,  0.1,  0.05,  0.01,  0.005,  0.001  and  0.0001  were  used.  For  this 
experiment,  the  following  NN  architecture  was  used:  540-input,  16-hidden  and  4-output  neurons. 
Phospholipid  fatty  acid  nomenclature.  Fatty  acids  are  designated  as  A:Ba)C,  where  A  is 
the  total  number  of  carbon  atoms,  B  is  the  number  of  double  bonds,  and  C  is  the  position  of  the 
double  bond  from  the  aliphatic  end  of  the  molecule.  Geometry  of  this  bond  is  indicated  'c'  for  cis 
and 't'  for  trans.  The  prefixes  'i'  and  'a'  refer  to  iso  and  anteiso  terminal  methyl-branching 
respectively  (26).  Mid-chain  methyl  branches  are  designated  by  'me'  preceded  by  the  position  of 
the  methyl  group  from  the  acid  end  of  the  molecule.  Cyclopropyl  fatty  acids  are  designated  as 
'ey'. 

Neural  network  software.  The  back-propagating  neural  network  program  used  in  this  study 
has  been  previously  described  (see  ref.  31).  A  binhex  and  binary  version  of  this  application  is 
available  through  an  anonymous  ftp  at  inlet.geol.sc.edu.  Additional  software  was  needed  to 
rapidly  calculate  the  hidden  layer  data  of  trained  NNs.  The  hidden  layer  data  were  used  for  cluster 
analysis.  A  stand-alone  application  for  the  Macintosh,  called  "hidden.app",  was  designed  by  using 
C++  software  (Symantec,  Release  5,  Cupertino,  CA).  A  binhex  and  binary  version  of  this 
application  is  also  available  at  our  ftp  site. 

Numerical  analysis.  To  evaluate  the  similarity  of  microbial  communities  from  sample  types 
and  worm  species,  cluster  analysis  was  performed,  using  the  Euclidean  distance  coefficient  and 
average  linkage  clustering  methods  (SPSS,  Incorp.).  The  robustness  of  group  memberships  was 
determined  by  training  the  NN,  calculating  the  hidden  layer  data,  and  performing  cluster  analyses. 
This  process  was  repeated  5  times,  resulting  in  the  production  of  5  dendrograms.  Comparisons  of 
cluster  membership  in  the  dendrogram  enabled  us  to  evaluate  robustness  of  group  memberships 
(15).  To  determine  which  of  the  PLFA  molecules  were  common  to  the  clusters,  the  mole 
percentage  of  PLFA  data  were  imported  into  Transform  3.01  software  (Spyglass,  Inc.,  Savoy, 

IL). 


RESULTS 

Optimization  of  the  neural  network.  The  effects  of  architecture  on  the  predictive  abilities  of 
NNs  were  determined  by  changing  the  number  of  hidden  layer  neurons  in  NNs  and  calculating  the 
percent  of  correct  predictions  for  both  sample  types  and  worm  species  (Table  1).  It  is  important  to 
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recognize  that  the  NNs  used  in  this  experiment  had  the  same  number  of  input  and  output  neurons 
but  different  numbers  of  hidden  layer  neurons.  Each  of  the  10  trained  NNs  used  to  evaluate 
predictive  power  of  NNs  were  trained  on  data  different  from  those  they  were  tested  on.  Moreover, 
each  cross-validation  procedure  was  performed  in  duplicate  to  provide  a  measure  of  the  prediction 
variability.  Increasing  the  number  of  hidden  layer  neurons  in  the  NN  decreased  the  number  of 
epochs  (i.e.,  number  of  completed  training  cycles)  required  to  reach  the  desired  error  tolerance 
(Table  1).  The  results  showed  that  the  number  of  neurons  in  the  hidden  layer  have  a  direct  effect 
on  the  learning  ability  of  NNs.  NNs  having  too  few  hidden  layer  neurons  (e.g.,  4  hidden  layer 
neurons)  did  not  converge  to  the  error  tolerance  of  0.01,  indicating  that  the  NN  was  not  able  to 
recognize  patterns  in  the  "fingerprints".  Although  increasing  the  number  of  neurons  in  the  hidden 
layer  decreased  the  mean  number  of  epochs  required  to  reach  the  desired  error  tolerance,  the  NN 
required  more  computation  time  to  reach  the  desired  error  tolerance  (data  not  shown). 

Increasing  the  number  of  hidden  layer  neurons  in  the  NNs  did  not  necessarily  improve  the 
number  of  correct  predictions.  As  indicated  in  Table  1,  eight  or  more  hidden  layer  neurons 
improved  the  number  of  correct  predictions  made  by  the  NNs.  However,  having  more  than  64 
hidden  layer  neurons  did  not  significantly  improve  NN  predictive  power  for  either  sample  type  or 
worm  species.  For  sample  type,  the  highest  number  of  correct  predictions  occurred  when  the 
number  of  hidden  layer  neurons  equalled  16  (Table  1).  For  worm  spjecies,  the  highest  number  of 
correct  predictions  occurred  when  64  hidden  layer  neurons  were  used  to  train  the  NN  (Table  1). 

For  subsequent  experiments,  the  following  architecture  was  used:  540-input,  16-hidden,  and  4- 
output  neurons. 

We  investigated  the  effects  of  different  error  tolerance  on  the  predictive  abilities  of  NNs  to 
establish  the  optimal  number  of  epochs  needed  for  training  NNs.  Decreasing  the  error  tolerance 
values  increased  the  number  of  epochs  required  to  train  the  NNs  (Fig.  lA)  and  therefore,  the 
amount  of  computational  time  (data  not  shown).  However,  the  percent  of  correct  predictions  did 
not  significantly  improve  beyond  an  error  tolerance  of  0.01  (Fig.  IB).  Presumably,  beyond  this 
threshold,  the  NNs  lost  their  ability  to  generalize  patterns  in  the  PLFA  profiles.  Based  on  these 
data,  the  error  tolerance  of  subsequent  experiments  was  set  at  0.01. 

Cross-validation  of  the  microbial  communities.  Examination  of  the  NN  predictions 
revealed  that  sample  type  (86.4  %)  was  more  often  correctly  predicted  than  worm  species  (63.6%), 
indicating  that  microbial  community  compositions  were  more  determined  by  whether  or  not  they 
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were  from  the  surface,  burrow  or  subsurface  than  a  particular  worm  species  (Table  2).  Incorrect 
predictions  provide  information  on  PLFA  profiles  having  common  characteristics.  The  overlap 
between  microbial  communities  predicted  for  burrow  and  subsurface  types,  and  burrow  and 
surface  types  are  presumably  due  to  aberrant  samples  (see  ref.  42)  since  there  was  little  or  no 
overlap  between  surface  and  subsurface  samples.  These  results  indicated  that  microbial 
communities  among  the  sample  types  are  distinctly  different  while  those  among  the  worm  species 
are  quite  similar  (Table  2). 

Cluster  analysis  of  the  microbial  community  PLFA  profiles.  A  dendrogram  showing 
the  relationship  between  microbial  communities  and  their  corresponding  PLFA  profiles  is  shown  in 
Figure  2.  The  dendrogram  represents  the  analysis  of  data  from  one  hidden  layer  of  one  NN  that 
was  trained  using  100%  of  the  PLFA  data.  We  compared  the  dendrogram  (Fig.  2)  to  those 
obtained  by  analyzing  the  data  from  four  additional  independently  trained  NNs  (data  not  shown). 
Although  weights  of  the  hidden  layer  changed  every  time  the  NNs  were  trained,  these  data  yielded 
similar,  but  not  identical  dendrograms,  indicating  that  group  memberships  for  most  samples  were 
robust.  Randomly-generated  samples  clustered  separately  from  the  rest  of  the  clusters  and  formed 
a  tight  cluster  at  a  rescaled  similarity  coefficient  of  ca.  87  in  all  dendrograms.  This  coefficient  was 
used  to  group  the  profiles  into  clusters. 

Based  on  the  composition  of  the  clusters  in  Figure  2,  group  memberships  appeared  to  be 
determined  by  both  sample  type  and  worm  species.  However,  the  other  dendrograms  support  the 
cross-validation  results  that  sample  type  has  a  much  greater  effect  on  cluster  membership  than 
worm  species.  Comparison  of  the  clusters  to  their  corresponding  PLFA  profiles  revealed  which 
PLFA  classes  and  which  individual  PLFAs  were  common  to  a  given  cluster  (Fig.  2)  and  allowed 
some  broad  conclusions  on  PLFA  distributions  to  be  reached.  Specifically,  surface  sediments 
characteristically  contained  high  levels  of  long  chain  (i.e.,  greater  than  or  equal  to  20  carbon  units 
in  length)  saturated  and  polyunsaturated  fatty  acids,  indicative  of  microeucaryotes  (14, 30, 33, 40, 
48, 49, 50, 52),  and  low  levels  of  terminally  and  mid-chain  branched  saturated  and  cyclopropyl 
fatty  acids,  indicative  of  low  levels  of  anaerobic  and/or  gram  positive  bacteria  (12, 13, 44, 47). 
Subsurface  sediments  showed  the  opposite  trend,  with  generally  low  levels  of  long  chain  fatty 
acids  and  high  levels  of  branched  fatty  acids.  Burrow  samples  had  PLFA  profiles  intermediate 
between  these  two  extremes,  reflecting  the  oxic/anoxic  interface  microenvironment  burrows  and 
tubes  represent.  Another  indication  of  this  is  the  relatively  high  abundance  of  cyclopropyl  fatty 
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acids,  particularly  in  B.  americam  tubes.  Cylcopropyl  fatty  acids  are  often  associated  with 
oxic/anoxic  interface  environments  (22). 


DISCUSSION 

In  this  study,  two  neural  computing  approaches  were  used  to  recognize  patterns  in  the 
PLFA  data  of  microbial  communities;  cross-validation  of  the  NNs'  predictions  and  cluster  analysis 
of  the  hidden  layer  of  trained  NNs.  To  our  knowledge,  this  is  the  first  study  to  use  two  different 
neural  computing  approaches  for  examining  complex  PLFA  patterns  of  natural  microbial 
communities.  We  employed  a  statistical  cross-validation  scheme  to  provide  information  on  the  true 
predictive  capabilities  of  NNs  and  cluster  analysis  of  the  hidden  layer  data  of  trained  NNs  to 
determine  the  relative  similarities  of  PLFA  profiles  among  the  samples.  This  study  and  a  previous 
study  (31)  have  shown  that  cluster  memberships  of  samples  were  robust  and  similar  regardless  of 
which  hidden  layer  was  used  for  cluster  analysis,  indicating  that  our  approach  was  valid.  The 
dendrograms  were  similar,  but  not  identical,  for  different  hidden-layer  data  because  the  weights 
change  every  time  the  NN  was  trained  and  the  weights  were  randomized  in  the  first  training  cycle. 
By  comparing  an  image  of  the  PLFA  data  to  its  corresponding  dendrogram  (Fig.  2),  we 
determined  which  patterns  in  the  microbial  communities  were  common  to  the  clusters,  providing 
insight  into  the  characteristics  likely  used  by  the  NN  to  classify  the  PLFA  data. 

This  study  also  used  the  cross-validation  approach  to  optimize  the  architecture  of  the  NN 
and  to  determine  the  extent  of  training  required  to  most  accurately  predict  the  sample  type  and 
worm  species.  The  NN  architecture  affected  the  mean  number  of  epochs  required  to  train  the  NN 
as  well  as  the  number  of  correct  predictions.  For  a  back-propagating  NN,  such  as  the  one  used  in 
this  study,  variation  in  the  NN  architecture  is  primarily  limited  to  changing  the  number  of  hidden 
layer  neurons  and/or  the  number  of  hidden  layers  since  the  number  of  input  and  output  neurons  are 
fixed  by  the  number  of  digits  in  the  raw  data.  We  used  one  hidden  layer  because  preliminary 
results  (this  study;  28)  showed  that  more  than  one  hidden  layer  was  rarely  beneficial.  A  low 
number  of  hidden  layer  neurons  increased  the  number  of  epochs  needed  to  train  the  NN  and 
decreased  the  percent  of  correct  predictions  presumably  because  the  NN  did  not  have  enough 
neurons  to  generalize  the  PLFA  profiles.  Similar  results  were  obtained  by  Weinstein  et  al.  (51). 
Having  more  than  64  hidden  layer  neurons  also  did  not  improve  the  percent  of  correct  predictions 
(Table  1),  indicating  that  the  NN  may  have  lost  its  ability  to  generalize.  Choosing  the  appropriate 
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number  of  hidden  layer  neurons  therefore,  is  extremely  important  for  optimizing  NN  performance. 


The  application  of  neural  computing  approaches  to  analyze  complex  data  is  relatively  new 
in  microbiology.  For  example,  NNs  have  been  used  to  identify  the  restriction  enzyme  patterns  of 
Escherichia  coli0\51'.\^  (10),  the  pyrolysis  mass  spectra  of  Mycobacterium  tuberculosis 
complex  species  (19),  the  promoter  sites  of  E.  coli  (23, 24),  protein-DNA  binding  sites  (32),  fatty 
acids  of  marine  heterotrophic  bacteria  (7),  and  stable  low-molecular  weight  rRNA  (31).  The 
advantage  of  using  neural  computing  approaches  is  that  they  allow  pattern  recognition  that  cannot 
be  substantiated  by  standard  statistical  methods  (28).  This  method  also  provides  numeric  and 
graphic  information  facilitating  interpretation  of  PLFA  data  for  the  analysis  of  natural  microbial 
communities. 
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FIGURE  LEGENDS 

Figure  1.  Relationship  of  error  tolerance  to  epochs  (A)  and  to  percent  correct  predictions  (B). 
Each  circle  represents  the  result  of  one  cross-validation  experiment.  Open  and  closed  circles  in  B 
represent  correct  predictions  for  sample  type  and  worm  species,  respectively. 

Figure  2.  A  dendrogram  showing  the  relationships  of  microbial  communities  from  collected 
samples  and  their  corresponding  PLFA  profiles.  The  similarity  scores  of  the  112  samples  were 
determined  by  calculating  the  squared  Euclidean  distance  and  clustering  the  data  using  average 
distance  linkage  methods.  The  23  randomly-generated  samples  are  represented  by  the  letter  R. 
Black,  gray  and  white  symbols  represent  burrow,  surface  and  subsurface  sample  types, 
respectively.  Circles,  squares  and  triangles  represent  Branchiasychus  americaria,  Balanoglossus 
aurantiacus,  and  Notomastus  lobatus,  respectively.  Each  PLFA  molecule  is  represented  by  a  box 
in  the  "fingerprint"  image,  whose  shade  is  based  on  phospholipid  fatty  acid  mol%  values;  white 
and  dark  boxes  representing  low  (0)  and  high  intensities  (>5),  respectively.  The  order  of  the 
PLFA  are  as  follows:  saturated,  14:0, 15:0,  16:0, 17:0, 18:0, 20:0, 22:0, 24:0;  branched  saturated, 
il4:0,  il5:0,  al5:0,  il6:0,  il7:0,  17: 1(o8c/al7:0;  mono  unsaturated  al6:0/16:l(D9c,  16:la)7c, 
16:la)7t,  16:lco5c,  16:la)13t,  17:la)6c,  18:lcD7t,  16:2(o6/brl5:0, 18:la)9c,  18:lco7c,  19:l(ol2c, 
19:1cd8c,  19:  1cd6c,  20:  l(o9c,  20:  1cd7c;  polyunsaturated.  Poly  17  (2+0)),  18:3a)6/10Mel7:0, 
18:4a)3/12Mel7:0,  18:2(d6,  18:3(o3/Brl7;l/il8:0, 20:4a)6,  20:5a)3,  20:3a)6,  20:4a)3,  20:3co3, 
22:5cd6,  22:6co3,  20:2a)6,22:4(o6,  22:5a)3,  poly  19  (3+(o),  Poly  22  (2(o); methyl-branched 
saturated,  10Mel6:0,  10Mel8:0, 1  lMel8:0;  cyclopropyl,  cy  17:0,  cyl9:0(a)7 ,8). 
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^  ABSTRACT 

2 

3  Infaunal  burrows  and  tubes  support  elevated  levels  of  microbial  biomass  and  activities 

4  relative  to  the  surrounding  bulk  sediments.  The  impacts  of  these  structures  on  sediment 

5  biogeochemistry  have  been  intensively  studied,  but  very  little  is  known  about  the 

6  microbiota  or  their  spatial  organization  within  these  microenvironments.  Distributions  of 

7  bacterial  cells  and  microcolonies  and  of  potentially  active  bacterial  cells  in  tubes  of  the 

8  onuphid  polychaete  Diopatra  cuprea  were  determined  using  confocal  scanning  laser 

9  microscopy.  Distributions  were  examined  among  different  tubes,  at  different  locations 

10  within  a  tube,  and  at  different  depths  in  the  biofilms  lining  these  tubes.  The  average 

1 1  quantity  of  cells  across  all  tubes  examined  was  5  .61  x  10*  per  cm^  and  there  was  no  trend 

12  in  the  distribution  of  cells  along  the  length  of  the  tube  or  with  increasing  depth  in  the 

13  biofilm.  Cell  size  across  all  tubes  collected  in  February  1997  averaged  0.085  pm^  but 

14  varied  seasonally.  Mean  cell  size  increased  with  increasing  depth  below  the  sediment 

15  surface  of  the  tube  site  sampled,  and  decreased  with  increasing  depth  in  the  biofilm. 

16  Microcolonies  were  abundant  and  not  confined  to  specific  depths  within  the  biofilm  or 

1 7  locations  within  the  tube.  Potentially  active  cells  were  observed  at  all  depths  within  the 

18  biofilm  and  at  least  46%  (mean  =  68%)  of  cells  at  any  given  site  were  potentially  active. 

19  We  hypothesize  that  the  observed  trends  in  cell  distribution  are  responsive  to  tube 

20  irrigation  by  the  host  animal. 


2 


INTRODUCTION 


1 

2 

3  Geochemical  cycling  of  the  major  biogenic  elements  in  nearshore  marine 

4  sediments  is  primarily  mediated  by  microorganisms.  However,  sediments  are  far  from 

5  homogeneous  and  some  specific  microenvironments  are  particularly  conducive  to 

6  microbial  activities.  Chemocline  environments,  such  as  those  found  at  the  sediment- 

7  water  interface,  are  well  known  to  be  sites  of  enhanced  microbial  activities  (Novitsky 

8  1983,  Novitsky  and  Karl  1986).  The  classic  model  of  these  interface  environments 

9  places  aerobic  metabolism  at  the  sediment  surface,  followed  by  various  anaerobic 

10  processes  in  deeper  sediment  layers  (e.g.  Aller  1988,  Fenchel  1969,  Revsbech  and 

1 1  Jorgensen  1986).  However,  surface  sediments  are  subject  to  high  levels  of  disturbance 

12  due  to  physical  and  biogenic  reworking  (Grant  1983,  Yingst  and  Rhoads  1 980). 

13  Sediment  mixing  can  greatly  stimulate  some  microbial  activities  by  increasing  the 

14  quantities  of  potential  electron  donors  and  acceptors  below  the  sediment- water  interface 

15  (Findlay  and  White  1983,  1984,  Findlay  et  al.  1990,  Krantzberg  1985).  Mixing  also 

16  disrupts  chemocline  microbial  communities  and  may  significantly  alter  the  interactions 

17  among  various  microbial  functional  groups  and  affect  their  potential  for  consortial 

18  metabolism  (Paerl  and  Pinckney  1996). 

19 

20  Macrofaunal  burrows  and  tubes  also  provide  significant  chemocline  niches  for 

21  microbial  growth.  These  structures  are  dug  into  anoxic  subsurface  sediments,  irrigated 

22  with  oxic  seawater,  and  provide  a  large  surface  area  for  diffusive  exchange  (Aller  1 980, 

23  Aller  and  Yingst  1978,  Boudreau  and  Marinelli  1994).  Irrigation  of  the  burrow  by  the 


3 


1  macrofaunal  host  organism  results  in  the  formation  of  radial  chemoclines  extending  out 

2  from  the  burrow  wall  into  the  surrounding  sediments  (Aller  1983, 1988,  Aller  and  Yingst 

3  1978,  Boudreau  and  Marinelli  1994,  Marinelli  and  Boudreau  1996).  These  radial 

4  chemoclines  provide  overlapping  gradients  of  potential  electron  donors  and  acceptors  and 

5  are  stable  over  time  spans  long  enough  for  development  of  complex  microbial 

6  commumties  (Aller  and  Aller  1 986,  Marinelli,  Wakeham,  Ringelberg,  White,  and  Lovell, 

7  unpublished).  Macrofaunal  burrows  vary  in  their  structural  properties,  ranging  from 

8  loosely  packed  sediment/mucopolysaccharide  aggregates,  as  seen  in  the  burrows  of  the 

9  terebellid  polychaete  Amphitrite  ornata,  to  well-defined  polysaccharide  and/or 

10  proteinaceous  tubes,  such  as  those  of  the  chaetopterid  polychaete  Chaetopterus 

1 1  variopedatus  and  the  onuphid  polychaete  Diopatra  cuprea  (Aller  1 983,  Aller  and  Yingst 

12  1978,  Myers  1972).  These  structures  promote  development  of  microbial  biofilms 

13  (Steward  et  al.  1996)  containing  high  levels  of  microbial  biomass  and  supporting  intense 

14  geochemical  cycling  activities  (Binnerup  et  al.  1992,  Kristensen  et  al.  1985, 1991). 

1 5  Previous  studies  have  shown  that  the  microbiota  of  burrow  lining  biofilms  differs  from 

1 6  that  of  adjacent  bulk  sediments  (Dobbs  and  Guckert  1 988,  Steward  et  al.  1 996)  and  that 

17  burrow  and  tube  microbiota  play  important  roles  in  the  mineralization  of  organic 

18  compounds  and  the  geochemical  cycling  of  biogenic  elements  (e.g.  Binnerup  et  al.  1992, 

1 9  Kristensen  et  al.  1 985,  1 99 1 ).  The  rates  of  geochemical  processes  within  the  burrow 

20  biofilm  rely  mainly  on  the  constituents  of  the  microbial  community  and  their  levels  of 

21  activity.  The  latter  are  obviously  highly  dependent  on  the  properties  of  the  burrow 

22  (stability  and  permeability)  and  on  the  irrigation  and  exudation  activities  of  the  host 

23  animal.  However,  the  distributions  and  potential  for  consortial  interactions  of  these 
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1  organisms  are  almost  completely  unknown  and  can  best  be  studied  using  microscopic 

2  analysis. 

3 

4  Recent  developments  in  microscopy,  specifically  confocal  scanning  laser 

5  microscopy,  allow  examination  of  bacterial  biofilms  in  a  non-invasive  manner  (Costerton 

6  etal.  1995).  Studies  to  date  have  primarily  examined  biofilms  grown  on  impermeable 

7  surfaces  and  little  is  known  about  natural  biofilms  growing  on  permeable  surfaces  and  the 

8  microbiota  inhabiting  them.  In  this  study  we  examined  distributions  of  total  and  active 

9  bacterial  cells  in  biofilms  lining  the  stable  tube  structure  constructed  by  a  ubiquitous 

10  infaunal  polychaete,  Diopatra  cuprea.  Several  key  features  of  the  tube  lining  biofilm 

1 1  microbiota  were  revealed. 

12 

1 3  CHARACTERISTICS  OF  THE  DIOPA  TEA  CUPREA  TUBE 

14 

15  The  D.  cuprea  tube  is  constructed  from  a  mucous,  sulphated  polysaccharide, 

16  released  from  secreting  glands,  which  hardens  when  it  comes  in  contact  with  seawater 

17  (Myers  1972).  Each  tube  displays  two  types  of  construction,  reinforced  and  unreinforced 

18  (Figure  1).  The  reinforced  section  rises  fi’om  the  sediment  surface  and  can  also  extend 

19  down  into  shallow  subsurface  sediments.  The  length  of  this  section  is  dependent  on  the 

20  degree  of  sediment  disturbance  and  the  animal  will  extend  it  if  the  exposed  portion 

21  becomes  buried.  The  reinforced  section  of  the  tube  is  decorated  with  debris  (small 

22  pebbles,  shell  fi-agments,  and  plant  detritus)  from  the  sediment  surface.  Occasionally, 
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1  tubes  will  have  two  reinforced  openings.  The  unreinforced  section  occurs  below  the 

2  sediment  surface  and  is  constructed  as  the  worm  moves  vertically  down  through  the 

3  sediments,  creating  a  mucous  sheath.  The  only  foreign  materials  that  may  be  attached  to 

4  this  portion  of  the  tube  are  sediment  particles  in  direct  contact  with  the  mucus  when  it  is 

5  secreted.  Additional  layers  of  mucus  are  laid  within  the  tube  over  time,  so  the  thickness 

6  of  the  tube  is  directly  proportional  to  its  age. 

7 

8  MATERIALS  AND  METHODS 

9 

10  Sampling  site.  Diopatra  cuprea  tubes  were  collected  from  an  intertidal  sandflat  located 

1 1  at  Oyster  Landing,  in  the  North  Inlet  saltmarsh  (79°W,  33°N),  near  Georgetown,  South 

12  Carolina,  USA.  Oyster  Landing  is  a  small  tidal  tributary  bordered  by 

13  alterniflora  marsh  and  oyster  beds  (Allen  et  al.  1992,  Lovell  and  Piceno  1994).  The 

14  sandflat  from  which  samples  were  collected  is  completely  exposed  during  low  tide  and 

15  consists  of  quartz  sand  (median  grain  size  of  376  pm)  containing  0.22%  silt  and  clay,  and 

16  0.83%  organic  matter,  by  dry  weight  (Lovell  and  Piceno  1994). 

17 

18  Sample  collection.  All  D.  cuprea  tubes  collected  for  this  study  were  built  by  solitary 

19  individuals  and  had  no  tube  branches.  The  amount  of  debris  and  firm  attachment  of  the 

20  debris  particles  to  the  tube  in  the  reinforced  section  made  it  impossible  to  moimt  samples 

21  of  intact  reinforced  segments  for  microscopic  examination  or  to  remove  the  debris 

22  without  disruption  of  the  inner  lining  biofilm.  Therefore,  only  the  unreinforced  sections 
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1  of  the  tubes  were  analyzed  in  this  study.  The  tube  caps  (apex  of  the  reinforced  section) 

2  were  cut  off  at  the  sediment  surface  and  the  tubes  exhumed  until  approximately  20-25  cm 

3  of  the  unreinforced  portions  were  exposed.  The  length  of  the  reinforced  portion  of  each 

4  tube  was  measured  and  the  approximate  distance  of  the  upper  end  of  the  unreinforced 

5  section  from  the  sediment  surface  determined.  From  each  tube  a  20  cm  sample  (the 

6  deepest  2  cm  of  the  reinforced  and  the  first  18  cm  of  the  unreinforced)  was  removed  and 

7  further  dissected  into  two  1 0  cm  segments,  designated  segment  A  and  B  (Figure  1 ). 

^  These  segments  were  selected  to  permit  examination  of  the  biofilm  at  close  proximity  to 
9  the  tube  inlet  and  at  greater  distance  from  it.  Segment  A  consisted  of  2  cm  of  the 

10  reinforced  and  the  first  8  cm  of  the  unreinforced.  Segment  B  was  the  remaining  1 0  cm  of 

11  the  unreinforced  tube  sample.  The  2  cm  reinforced  end  of  segment  A  and  the  first  2  cm 

1 2  of  segment  B  permitted  convenient  handling  of  the  tube  samples  without  disruption  of  the 

13  biofilm  in  the  remaining  8  cm  segments. 

14 

15  Five  D.  cuprea  tubes  were  collected  in  February  1 997.  The  water  and  sediment 

16  temperatures  at  the  time  of  collection  were  20°  C  and  19°  C,  respectively.  The  samples 

17  were  placed  in  0.2  pm  filtered  2.5%  fonnalin,  transported  to  Columbia,  SC,  and  stored  at 

18  4°  C  in  the  dark  for  no  more  than  one  week  prior  to  analysis.  Five  additional  D.  cuprea 

19  tubes  were  collected  in  October  1997.  The  water  temperature  was  26°  C  and  the 

20  sediment  temperature  was  25°  C  at  that  time.  Tubes  were  handled  as  described  above, 

21  but  with  the  following  modifications:  (1)  only  segment  A  for  each  tube  was  collected 

22  and  (2)  tube  samples  were  not  formalin  fixed.  These  samples  were  transported  to 

23  Columbia,  SC  and  held  in  a  recirculating  sea  water  tank  overnight. 
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1  Sample  preparation  and  staining  procedures.  Each  ofthe  ten  February  1997  tube 

2  segments  was  secured  onto  a  glass  microscope  slide  with  rubber  bands  approximately  5 

3  cmapmt.  The  2  cm  ends  designated  for  handling  were  then  removed.  The  portion  of  the 

4  8  cm  segment  between  the  two  rubber  bands  was  cut  longitudinally  and  unrolled  with  the 

5  inner  tube  surface  facing  up.  These  samples  were  stained  for  eniuneration  and 

6  measurement  of  total  cells  and  microcolonies. 

7 

8  SYTOX  Green  was  selected  as  the  primary  fluorescent  stain  after  substantial 

9  background  fluorescence  problems  were  encountered  with  acridine  orange  and 

10  fluorescein  isothiocyanate.  DAPI  (4,6-diamidino-2-phenylindole)  was  not  tested  since 

11  the  excitation  wavelength  needed  for  this  dye  could  not  be  obtained  with  our  confocal 

12  microscope  system.  Fluorescence  background  from  SYTOX  Green,  SYTO-9,  and 

13  propidium  iodide  was  quite  low  and  permitted  detailed  assessment  of  cell  distributions 

14  and  appearances.  One  hundred  pL  of  a  12.5  pM  solution  of  SYTOX  Green  (Molecular 

15  Probes,  Inc.,  Eugene,  OR)  were  added  directly  onto  the  inner  surface  and  the  sample 

16  incubated  at  room  temperature  in  the  dark  for  5  minutes.  SYTOX  Green  penetrates  only 

17  fixed  cells  (Roth  et  al.  1997),  complexing  with  DNA.  The  stain-DNA  complex 

18  fluoresces  green  when  excited  (excitation/emission;  504  nm/523  nm).  After  incubation,  a 

1 9  solution  of  1 ,4-diazobicyclo[2,2,2]octane:glycerol  ( 1 :3)  (Sigma,  St.  Louis,  MO)  was 

20  added  to  reduce  photobleaching.  A  22  mm  x  40  mm  coverslip  was  taped  down  on  top  of 

21  the  sample,  holding  it  flat.  The  unfixed  October  1997  tube  segments  were  mounted  as 

22  described  above  and  stained  for  the  enumeration  of  active  and  inactive  cells  (Lawrence  et 

23  al.  1 997).  Metabolically  active  bacteria  create  a  proton  gradient  across  their  plasma 
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1  membrane  through  active  transport  of  protons  out  of  the  cell.  This  electrochemical 

2  gradient,  called  the  proton  motive  force,  or  membrane  potential,  provides  a  source  of 

3  energy  that  can  be  used  to  drive  ATP  synthesis.  SYTO-9  (excitation/emission:  488 

4  nm/509  nm;  Molecular  Probes,  Inc.,  Eugene,  OR)  only  penetrates  cells  having  membrane 

5  potential  (Lebaron  et  al.  1998)  and  fluoresces  green  upon  excitation.  These  cells  are 

6  considered  to  be  potentially  active.  Propidium  iodide  (excitation/emission:  535  nni/617 

7  nm;  Molecular  Probes,  Inc.,  Eugene,  OR)  only  penetrates  cells  lacking  membrane 

8  potential  (Kaneshiro  et  al.  1993)  and  fluoresces  red  upon  excitation.  These  cells  are 

9  considered  to  be  inactive.  One  hundred  pL  of  a  50  pM  SYTO-9, 149  pM  propidium 

10  iodide  staining  solution  were  added  to  the  samples,  which  were  incubated  and  mounted  as 

1 1  described  above. 

12 

13  Confocal  scanning  laser  microscopy.  A  Bio-Rad  MRC  -  1000  Laser  Scanning 

14  Confocal  Microscope  (Hercules,  CA)  equipped  with  a  krypton/argon  ion  laser  was  used 

15  to  produce  fluorescence  images  of  the  tube  lining  biofilms.  Confocal  microscopy  allows 

16  optical  sectioning  of  the  biofilm  using  an  adjustable  z-axis  stepping  motor.  At  set 

17  successive  intervals,  the  tunable  laser  (set  to  the  appropriate  excitation  wavelength)  is 

18  scanned  over  the  sample  and  excites  fluorescent  dye  complexes  in  a  specific  focal  plane, 

19  allowing  detection  of  stained  cells  at  that  specific  depth  (z)  in  the  biofilm.  Data  were 

20  collected  at  1  pm  depth  intervals  in  the  biofilm  until  the  tube  wall  was  encountered, 

21  providing  a  complete  z-profile  of  images  from  the  top  of  the  biofilm  to  the  tube  wall. 

22  The  laser  does  not  penetrate  below  the  tube  wall  inner  surface,  and  no  fluorescence  is 

23  observable  below  that  depth  (Figure  1). 
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1 

2  Two  z-profiles  were  taken  at  each  of  five  sites  in  each  of  the  SYTOX  Green 

3  stained  tube  segments.  One  profile  of  each  pair  was  collected  at  1 OOOX  magnification  to 

4  allow  determination  of  single  cell  distribution,  the  other  at  400X  magnification  for 

5  determination  of  microcolony  distribution.  The  first  site  in  each  tube  segment  was 

6  approximately  2.5  cm  fi-om  the  top  of  the  segment  (site  1),  with  successive  sites  (sites  2- 

7  5)  at  5  mm  intervals  along  the  length  of  the  segment.  The  dimensions  of  the  fields  of 

8  view  were  dependent  on  the  objective  used  (lOOOX:  1 14  x  76  pm,  400X:  288  x  192  pm). 

9  Although  data  sets  for  single  cell  and  microcolony  distributions  were  collected  at 

10  approximately  the  same  sites,  each  is  considered  to  be  an  independent  sample  due  to 

1 1  microscale  differences  in  biofilm  surface  texture  and  the  difference  in  the  areas  analyzed 

12  using  lOOOX  and  400X  magification. 

13 

14  Profiles  were  also  collected  at  five  sites  at  lOOOX  magnification  in  each  of  the 

1 5  activity  stained  samples.  A  split  spectral  line  (488  nm  and  568  nm)  was  used  to  excite 

16  both  SYTO-9  and  propidium  iodide  in  order  to  visualize  the  potentially  active  (green) 

1 7  and  non-active  (red)  cells  simultaneously.  Data  were  not  collected  on  microcolonies  in 

18  segments  stained  to  determine  potentially  active  bacteria  due  to  the  bright  fluorescence  of 

19  SYTO-9,  which  tended  to  overwhelm  the  fluorescence  of  propidium  iodide  in 

20  multicellular  aggregates. 

21 

22  Data  collection.  Data  collected  at  1  pm  depth  intervals  in  the  biofilm  to  determine  the 

23  distribution  of  single  bacterial  cells  (February  1997  tubes)  were  compiled  into  a  series  of 
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1  5  |im  intervals  (i.e.  0-5  ^m,  6-10  ^im,  11-15  |j.m,  etc.)  and  the  number  of  cells  present  in 

2  each  interval  counted.  A  total  of  50  sites  were  examined  for  single  cell  distribution 

3  within  the  five  tubes,  but  data  for  one  site  in  one  tube  (tube  1,  segment  A,  site  1)  were 

4  lost  due  to  faulty  digital  storage.  In  order  to  avoid  duplicate  counting  of  cells  that  were  at 

5  the  border  of  two  adjacent  5  pm  intervals,  each  interval  image  (except  0-5  pm)  of  a 

6  sequential  pair  of  intervals  was  highlighted  in  a  different  color  (red  or  green)  and  the  two 

7  images  were  overlapped.  Cells  present  in  both  intervals  were  indicated  by  yellow 

8  highlighting  and  were  counted  only  in  the  first  of  the  two  intervals.  Bacterial  cells 

9  present  in  multicellular  aggregates  of  more  than  5  cells  were  not  counted  in  these 

10  profiles.  Lengths  and  widths  of  10  arbitrarily  selected  cells  were  determined  for  each 

1 1  depth  interval  by  optical  magnification  of  collected  images  and  measurement  using 

12  calipers.  If  less  than  10  cells  were  present  in  an  interval,  then  each  cell  was  measured. 

13  Mean  cell  size  (average  biovolume  of  a  single  cell)  for  each  interval  was  calculated  from 

14  the  length  and  width  measurements  determined  for  single  cells  (Norland  1993).  The  total 

15  biovolume  of  single  cells  for  each  5  pm  interval  was  calculated  by  multiplying  the 

16  number  of  cells  in  that  interval  by  the  mean  cell  size  for  that  interval. 

17 

1 8  Fifty  sites  were  analyzed  for  the  distribution  of  all  types  of  multicellular 

19  aggregates.  All  multicellular  aggregates,  whether  tightly  packed  and  representing 

20  relatively  high  densities  of  cells,  or  loosely  arranged,  with  relatively  few  cells,  were 

21  designated  as  microcolonies.  For  each  site,  the  1  pm  images  collected  at  400X 

22  magnification  were  compiled  in  3  pm  intervals.  Some  sites  did  not  contain 

23  microcolonies.  In  some  others,  profiles  were  not  completely  perpendicular  to  the  tube 
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1  wall,  prohibiting  accurate  measurement  of  depths  in  the  biofilm.  A  total  of  46  sites 

2  yielded  useful  data  on  microcolony  distributions  and  sizes.  Microcolony  cross  sectional 

3  area  in  each  of  the  3  pm  intervals  was  measured  using  MorphoSys  (Meacham  and 

4  Duncan  1991).  Biovolumes  of  the  microcolonies  were  estimated  by  multiplying  the 

5  microcolony  cross  sectional  areas  in  the  sequential  3  pm  intervals  by  the  total  depth  of 

6  biofilm  through  which  the  colonies  extended.  No  attempt  was  made  to  differentiate 

7  biovolume  occupied  by  cells  from  that  containing  only  exopolymeric  materials.  Due  to 

8  microscale  differences  in  the  sites  at  which  single  cell  and  microcolony  data  were 

9  collected,  total  biovolume  of  bacteria,  which  would  include  eill  cells  regardless  of 

10  differing  modes  of  growth,  were  not  calculated. 

11 

12  Depth  intervals  (5  pm)  of  each  profile  at  the  5  sites  in  each  of  the  tube  segments 

13  stained  to  determine  potentially  active  bacteria  were  compiled  as  described  above  and  the 

14  number  of  potentially  active  cells  and  non- viable  or  inactive  cells  in  each  of  the  resulting 

15  5  pm  depth  intervals  were  coimted.  Duplicate  counting  of  cells  located  at  the  border  of 

16  two  adjacent  intervals  was  avoided  by  aligning  the  two  intervals  side  by  side  and 

17  comparing  cell  placement  in  the  field  of  view.  One  arbitrarily  selected  profile  from  each 

18  activity  stained  tube  segment  was  used  for  determination  of  average  sizes  of  active  and 

19  inactive  cells  (tube  1 ,  site  1 ;  tube  2,  site  1 ;  tube  3,  site  5;  tube  4,  site  2;  tube  5,  site  3)  as 

20  described  above.  Lengths  and  widths  of  all  potentially  active  and  inactive  cells  were 

21  measured  and  mean  cell  sizes  calculated  for  each  of  the  5  pm  intervals  in  these  profiles  as 

22  described  above. 

23 
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1  Statistical  analyses:  The  SAS  System  for  Windows  (version  6.12,  SAS  Institute,  Inc., 

2  Cary,  NC,  USA)  was  used  to  perform  all  statistical  analyses.  Nested  Analysis  of 

3  Variance  was  used  to  determine  significant  differences  in  cell  numbers,  cell  size,  and 

4  biovolume  among  tubes,  among  segments,  and  among  sites  for  data  collected  from  D. 

5  cuprea  tubes  stained  with  SYTOX  Green.  Regression  analysis  was  used  to  evaluate 

6  relationships  between  cell  numbers,  mean  cell  size,  or  biovolume,  and  depth  in  the 

7  biofilm  for  all  sets  of  tube  segments.  Cell  size  variance  was  modeled  using  the  gamma 

8  distribution  with  Generalized  Linear  Model  (GLIM),  using  Generalized  Estimating 

9  Equations  (GEE)  to  accommodate  the  nested  design.  All  data  were  tested  for  normality 

10  and  statistical  tests  were  performed  using  a  significance  level  of  a  =  0.05. 

11 

12  RESULTS  AND  DISCUSSION 

13 

14  Bacterial  single  cell  distributions  and  sizes 

15 

16  Confocal  microscopy  revealed  that  the  tubes  of  Diopatra  cuprea  supported 

17  substantial  biofilms  containing  abimdant  bacteria  occurring  as  single  cells  eind  in 

18  microcolonial  formations.  The  thickness  of  D.  cuprea  tube  biofilm  ranged  from  30  pm  to 

19  110  pm,  with  an  average  thickness  of  60  pm.  The  quantity  of  individual  bacterial  cells  at 

20  a  given  site  ranged  from  1.38  x  10*  to  1.33  x  10^  cells*cm‘^  (mean  =  5.61  x  10* 

21  cells«cm*^.  Figure  2 A).  Nested  analysis  of  variance  revealed  that  cell  numbers  differed 

22  significantly  from  site  to  site  (df  8,  3 1 ;  F  =  3.00;  p  <  0.02),  but  not  among  tubes  (df  4,4;  F 

23  =  3.14;  p>  0.14)  or  between  segments  A  and  B  (df  4,8;  F  =  1.01;  p>  0.45).  There  was 

24  no  consistent  pattern  of  cell  distribution  >vithin  the  tube  lining  biofilm  (Figure  3). 
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1  Twenty-five  of  the  49  sites  had  highest  cell  numbers  within  1 0  ^im  of  the  biofilm  surface 

2  Twenty-three  had  additional  subsurface  peaks  in  cell  numbers.  Cell  numbers  decreased 

3  significantly  (regression  analysis,  p  <  0.05  in  all  cases)  with  increasing  depth  in  the 

4  biofilm  at  28  of  the  49  sites,  but  the  remaining  sites  showed  no  significant  changes  with 

5  depth. 

6 

7  Unlike  cell  numbers,  mean  bacterial  cell  sizes  did  change  with  depth  below  the 

8  sediment  surface.  Cell  sizes  at  different  sites  in  the  tube  segments  ranged  from  0.032 

9  pm^  to  0.285  pm^  (grand  average:  0.085  pm^.  Figure  2B).  Nested  analysis  of  variance 

10  revealed  significant  differences  in  mean  cell  size  between  segments  A  and  B  (df  4,8;  F  = 

1 1  4.93;  p  <  0.03),  with  cells  in  segment  B  averaging  larger  than  those  in  segment  A.  No 

12  significant  differences  in  mean  cell  size  were  found  among  tubes  (df  4,4;  F  =  0.60;  p  > 

1 3  0.68)  or  sites  (df  8,31;F=1.19;p>  0.3).  Based  on  D.  cuprea  behavior  and  bacterial 

14  stress  responses,  this  trend  is  quite  different  fi-om  what  would  be  expected.  D.  cuprea 

1 5  spends  most  of  its  time  in  the  upper  portion  of  the  tube  near  the  sediment  surface  and 

16  irrigation  of  the  tube  is  frequent  but  not  continuous.  At  the  depths  of  the  tube  sites 

17  examined  (10  cm  to  36  cm  below  the  sediment  surface),  the  only  oxygen  available  to  the 

18  microbiota  and  a  substantial  portion  of  the  nutrient  inputs  are  provided  by  host  irrigation. 

19  Consequently,  sites  at  greater  depths  should  experience  greater  oxygen  and/or  nutrient 

20  limitation  due  to  the  lesser  impact  of  irrigation.  Bacteria  in  nature  that  are  severely 

21  nutrient  limited  are  typically  smaller  than  their  more  resource-sufficient  coimterparts 

22  (Chesbro  et  al.  1990,  Kaprelyants  et  al.  1993,  Kjelleberg  et  al.  1983).  If  cells  located 
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1  deeper  in  the  tube  were  oxygen  and/or  nutrient  limited  they  should  have  been  smaller 

2  than  those  at  shallower  depths,  but  that  was  not  the  case. 

3 

4  There  was  also  a  clear  pattern  in  the  distribution  of  bacterial  cell  size  with  depth 

5  in  the  burrow  lining  biofilm  (Figure  3).  The  largest  mean  cell  sizes  were  found  in 

6  intervals  within  10  pm  of  the  biofilm  surface  at  32  of  the  49  sites.  Nineteen  sites  had 

7  subsurface  peaks  in  cell  size.  Almost  all  sites  (42  out  of  49)  showed  significant 

8  decreases  in  mean  cell  size  (regression  analysis,  p  <  0.05  in  all  cases)  and  of  variance  in 

9  cell  size  (E  =  -0.0005;  SE  =  0.0002;  Z  =  -1 .906;  p  <  0.01)  with  increasing  depth  into  the 

10  biofilm.  This  trend  may  have  been  due  to  nutrient  or  oxygen  limitation  at  the  base  of  the 
1  ]  biofilm. 

12 

13  The  biovolume  of  single  cells  at  a  given  site,  a  fimction  of  the  grand  average  of 

14  cell  sizes  (i.e.  average  of  mean  cell  sizes  for  all  intervals  at  that  site)  and  cell  numbers, 

15  ranged  from  7.26  x  10^  to  3.53  x  10*  pm*  •  cm'*  (mean  =  5.33  x  lO’  pm*  •  cm'*.  Figure 

16  2C).  Nested  analysis  of  variance  showed  that  biovolume  differed  significantly  among 

17  sites  (df  8,3 1 ;  F  =  2.39;  p  <  0.04)  but  not  between  segments  A  and  B  (df  4,8;  F  =  2.06;  p 

18  >  0. 1 7)  or  among  tubes  (df  4,4;  F  =  0.95;  p  >  0.5 1 ).  The  highest  biovolumes  in  24  of  49 

19  sites  were  within  10  pm  of  the  biofilm  surface  while  12  sites  had  subsurface  bio  volume 

20  peaks.  There  was  a  significant  relationship  between  biovolume  and  depth  in  the  biofilm 

21  at  34  of  49  sites,  and  bio  volume  decreased  with  increasing  depth  in  the  biofilm  in  at  least 

22  half  of  the  sites  in  each  tube  (Figure  2). 

23 
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1  Microcolony  distribution 

2 

3  Microcolonies  covered  a  broad  range  of  sizes  and  were  found  throughout  the 

4  length  of  the  D.  cuprea  tube  and  at  all  depths  within  the  biofilm.  The  number  of 

5  microcolonies  at  a  given  site  ranged  fi-om  1  to  3 1  (mean  =  5)  and  these  formations  were 

6  located  at  various  depths  throughout  the  biofilm.  Microcolony  biovolume  ranged  from 

7  16  fj,m^  to  2.5  X  10^  pm^  (mean  =  1 .28  x  10"*  pm^).  Nested  analysis  of  variance  revealed 

8  that  the  microcolony  biovolumes  did  not  differ  significantly  among  sites  (df  8,28;  F  = 

9  1.55;  p  >  0.18),  between  segments  A  and  B  (df  4,8;  F  =  1.83;  p  >  0.21),  or  among  tubes 

10  (df  4,4;  F  =  1 .96;  p  >  0.26).  The  structures  of  microcolonies  varied  from  very  tightly 

1 1  organized  formations  to  loose  aggregations  of  cells  and  expolymers.  The  potential  of 

12  these  structures  to  facilitate  interactions  among  different  bacterial  species  (Paerl  and 

13  Pinckney  1996)  is  clear  and  their  formation  may  be  very  important  to  the  overall 

14  activities  of  the  biofilm. 

15 

16  Potentially  active  vs.  inactive  cells 

17 

18  D.  cuprea  tubes  were  also  stained  to  permit  determination  of  the  proportions  of 

19  cells  that  were  intact  and  potentially  active  versus  those  lacking  membrane  potential  and 

20  thus  presumably  inactive  (Lawrence  et  al.  1997).  The  total  number  of  cells  per  site  in 

21  segments  stained  to  determine  bacterial  activity  ranged  fi'om  2.62  x  1 0*  to  1 .71  x  1 0’  cells 

22  •cm  (mean  =  6.1 1  x  10  •cm' ).  The  proportion  of  total  cells  that  were  potentially 

23  active  was  quite  high  (Figure  4),  ranging  from  0.46  to  0.85  (mean  =  0.68),  and  did  not 
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1  significantly  differ  among  sites  (df  4,16;  F  =  1.33;  p  >  0.29).  Potentially  active  cells  were 

2  observed  at  all  depths  within  the  biofilm  (Figure  5)  and  the  proportion  of  cells  that  were 

3  potentially  active  was  similar  at  all  depths  at  19  of  the  25  sites.  The  proportion  of  active 

4  cells  increased  with  increasing  depth  in  the  biofilm  at  three  sites  and  decreased  with 

5  increasing  depth  at  the  remaining  three  sites  (p  <  0.05  in  all  cases).  The  abundance  of 

6  potentially  active  cells  was  not  surprising  since  previous  studies  have  shown  that  infaimal 

7  burrows  are  sites  of  intense  microbial  activity  and  contain  substantial  amounts  of  viable 

8  microbial  biomass  (Aller  1980,  Aller  and  Yingst,  1978,  Dobbs  and  Guckert  1988, 

9  Steward  et  al.  1996).  What  was  surprising  was  that  cells  at  depth  in  the  biofilm  and  at 

10  sites  far  below  the  sediment  surface  appear  to  be  as  capable  of  activity  as  those  near  the 

1 1  biofilm  surface  and  at  sites  near  the  sediment  surface. 

12 

13  Sizes  of  potentially  active  and  inactive  cells  were  also  measured  at  one  site  in 

14  each  of  the  five  tube  segments.  The  average  size  of  potentially  active  cells  ranged  from 

15  0.19  jim^  to  0.42  pm^  (grand  mean  =  0.3  pm^)  and  was  significantly  larger  than  that  of 

16  inactive  cells  (range:  0.15  pm^  to  0.28  pm^,  grand  mean  =  0.18  pm^)  across  all  five  tubes 

17  (df  4;  T  =  4.37;  p  <  0.01)  (Figure  6).  There  was  no  general  pattern  in  the  distribution  of 

18  potentially  active  or  inactive  mean  cell  sizes  in  relation  to  the  depth  in  the  biofilm  (Figure 

19  7).  The  variance  of  potentially  active  cell  size  significantly  decreased  with  increasing 

20  depth  in  the  biofilm  (E  =  -0.00 1 7;  SE  =  0.0006;  Z  =  -2.945 ;  p  <  0.0 1 )  while  the  variance 

21  of  inactive  cell  size  increased  with  depth  (E  =  0.0005;  SE  =  0.0004;  Z  =  -1 .416;  p  <  0.01). 

22  The  observed  size  differences  between  potentially  active  and  inactive  cells  were 
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1  consistent  with  reports  of  stressed  or  starved  cells  averaging  smaller  than  actively 

2  growing  cells  of  the  same  species  (reviewed  in  Kjelleberg  et  al.  1987). 

3 

4  Comparisons  were  also  made  between  mean  cell  numbers  and  sizes  in  the  A 

5  segments  of  tubes  collected  in  October  1997  and  February  1997.  There  was  no 

6  significant  difference  in  the  mean  number  of  cells  per  tube  between  the  two  sets  (df  1 ,8;  F 

7  =  0.3 1 ;  p  >  0.50).  However,  there  was  a  difference  in  mean  cell  size.  Mean  cell  sizes  in 

8  tubes  collected  in  October  1997  for  activity  staining  were  significantly  larger,  0. 1 8  -  0.38 

9  pm^,  than  those  in  tubes  collected  in  February  1997, 0.058  -  0.09  pm^  (df  1,8;  F  =  0.093; 

10  p<0.01).  One  clear  difference  between  the  February  and  October  samples  was 

1 1  formaldehyde  fixation  of  the  February  samples.  Turley  and  Hughes  (1 992)  reported 

12  about  10%  reduction  in  cell  size  (estimated  from  digitized  cell  images)  in  fixed  samples 

13  stored  for  39  to  53  days,  but  provided  no  data  for  shorter  storage  times.  Troussellier  et  al. 

14  (1995),  using  flow  cytometric  light  scatter  techniques,  found  no  measurable  changes  in 

15  cell  structure  or  morphology  in  samples  fixed  for  a  much  shorter  time  period  (24  h).  The 

16  short  times  (<1  week)  our  fixed  samples  were  stored  prior  to  analysis  seem  unlikely  to 

17  produce  significant  cell  shrinkage.  Even  shrinkage  comparable  to  that  reported  by  Turley 

18  and  Hughes  (1992)  could  not  account  for  the  almost  three-fold  difference  in  mean  cell 

19  size  we  observed  between  the  February  and  October  samples.  A  more  likely  explanation 

20  for  the  difference  in  cell  size  can  be  based  on  the  biology  of  D.  cuprea.  The  water 

21  temperature  in  October  was  warmer  than  in  February  and  the  polychaetes  likely  irrigated 

22  more  actively  in  October  (e.g.  Mangum  and  Sassaman  1969).  More  frequent  irrigation 

23  during  the  warmer  period  would  increase  the  availability  of  oxygen  and  nutrients. 


18 


1  Warmer  water  temperature  would  also  have  a  positive  effect  on  the  bacterial 

2  communities,  resulting  in  faster  grotvth  and  an  increase  in  mean  cell  size. 

3 

4  The  differences  in  the  relationship  of  mean  cell  size  and  depth  in  the  biofilm 

5  between  the  two  collection  dates  could  also  be  a  result  of  irrigation  intensity  and 

6  frequency.  In  February  cells  became  smaller  and  more  uniform  in  size  with  increasing 

7  depth  in  the  biofilm.  This  trend  was  not  observed  in  tubes  collected  in  October  and  could 

8  be  due  to  oxygen  and  nutrient  limitation  of  microorganisms  located  deeper  in  the  biofilm. 

9 
10 

^  ^  Recent  studies  of  many  laboratory-cultivated  biofilms  have  revealed  channels  of 

12  low-density-exopolymer  throughout  the  material  (e.g.  Stewart  et  al.  1993,  Stoodley  et  al. 

13  1994,  Wolfaardt  et  al.  1 994).  These  channels  are  thought  to  be  essential  for  nutrient 

14  transport  to  organisms  deep  within  the  biofilm  (Costerton  et  al.  1 995,  Lawrence  et  al. 

1 5  1 994)  and  can  be  revealed  through  staining  with  fluorescently  tagged  lectins,  which  bind 

16  specifically  to  the  biofilm  polysaccharides.  Staining  with  fluorescein  isothiocyanate- 

1 7  labeled  concanavalin  A  (Sigma,  St.  Louis,  MO),  which  selectively  binds  a- 

1 8  mannopyranosyl  and  a-glucopyranosyl  residues  in  polysaccharides,  revealed  no  channels 

1 9  within  the  D.  cuprea  tube  biofilm  (data  not  shown).  D.  cuprea  deposits  many  layers  of 

20  mucus  over  time  and  cells  that  are  on  the  surface  of  the  tube  lining  eventually  become 

21  embedded  within  the  tube  matrix.  Any  channels  that  might  arise  could  be  filled  with 

22  worm  mucus  on  the  next  cycle  of  polymer  deposition.  However,  channels  may  not  be 

23  required  for  active  microbial  metabolism  throughout  the  tube  lining  biofilm.  Aller  (1983) 

24  suggested  that  the  tube  linings  of  macro-infaunal  organisms  could  be  considered  to  be 

25  highly  hydrated  porous  membranes  (Aller  1 983).  Diffusion  of  small  molecules  through 
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1  the  tube  (including  the  biofilm)  would  not  be  hindered  by  such  materials  and  thus  activity 

2  of  bacteria  at  depth  in  the  biofilm  could  be  supported. 

3 

4  The  D.  cuprea  tube  supports  an  abundant  and  varied  microbial  community  that 

5  contains  an  unusually  high  proportion  of  potentially  active  cells.  These  findings  are 

6  consistent  with  the  high  levels  of  microbial  biomass  and  activity  reported  for  various 

7  infaunal  burrows  and  tubes.  The  trends  observed  in  cell  size  distributions  are  very 

8  intriguing  in  that  they  may  reflect  the  nutritional  and/or  stress  regimes  the  organisms 

9  experience  in  response  to  irrigation  of  the  tube  by  D.  cuprea.  Manipulative  experiments 

1 0  exploring  the  impact  of  tube  irrigation  on  biofilm  microbiota  and  molecular  biological 

1 1  characterization  of  these  organisms  are  ongoing. 

12 
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Figure  Legends 

.  Ci 

2 

3  Figare  1 .  Schematic  diagram  of  Diopatra  aiprea  tube  stmcture  and  the  sampling 

4  procedure  used  to  deteimine  distributioiis  of  bacteria  in  D.  cuprea  tube  lining  biofilms. 

5  For  each  tube,  two  8  cm  segments  of  the  unreinforced  portion  (segment  ^  and  B)  were 

6  collected.  In  each  of  the  segments,  five  sites  (1-5)  were  examined.  At  each  site,  confocal 

7  images  were  collected  at  1  pm  depth  intervals  fi-om  the  surface  of  the  biofilm  to  the  tube 

8  wall. 

9 

10  Figure  2.  Cells  per  cm’  (A),  mean  cell  size  (B).  and  biovolumc  (cell  numbers  multiplied 

11  by  mean  cell  size  at  each  site)  (C)  in  DIopaIra  cuprea  tube  lining  biofilms.  Mean  (±  SD) 

12  is  shown  for  each  of  the  five  sites  within  each  segment  across  all  5  tubes  (n  =  5  except  for 

1 3  segment  A,  site  1 ,  where  n  =  4). 

14 

15  Figure  3.  Distribution  of  cell  numbers  (bars),  mean  cell  size  (±  SD,  n  =  1-10)  (closed 

16  circles)  and  biovolume  (cell  numbers  multiplied  by  mean  cell  size  per  5  pm  interval) 

17  (open  triangles)  with  increasing  depth  into  Diopatra  cuprea  mbe  lining  biofihn  for  sites  1 

1  *  (A),  2  (B),  3  (C),  4  (D)  and  5  (E)  in  segment  A  of  tube  2, 

19 

20  Figure  4.  Total  (black  bars)  and  potentially  active  (open  bats)  cells  per  cm’  for  each  site 

21  in  segment  A  of  Diopatra  cuprea  tube  lining  biofilms.  Mean  (+  SD)  is  shown  for  each 

22  site  for  5  tubes. 

23 
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1  Figure  5.  Total  (black  bars)  and  potentially  active  (open  bars)  cells  per  cm^  in  segment 

2  A  of  Diopatra  cuprea  tube  lining  biofilms.  Mean  (±  SD)  at  each  depth  in  the  biofilm, 

3  across  all  sites,  is  shown.  The  number  of  samples  varies  from  1  to  25  at  different  depths 

4  in  the  biofilm  due  to  differences  in  biofilm  thickness  at  different  sites. 

5 

6  Figure  6.  Potentially  active  (closed  circles)  and  inactive  (open  triangles)  cell 

7  sizes  in  Diopatra  cuprea  tube  lining  biofilms.  Mean  (±  SD)  is  shown  for  one  site  in  each 

8  tube.  Sample  size  (n)  for  potentially  active  cells  ranges  from  43  to  94;  sample  size  (n)  for 

9  inactive  cells  ranges  from  33  to  79. 

10 

1 1  Figure  7.  Mean  potentially  active  (closed  circles)  and  inactive  (open  triangles)  cell  sizes 

12  in  Diopatra  cuprea  tube  lining  biofilms.  Means  (±  SD)  at  each  depth  in  the  biofilm, 

1 3  across  all  five  sites,  is  shown.  The  number  of  samples  at  a  given  depth  in  the  biofilm 

14  varies  due  to  differences  in  biofilm  thickness  (n  =  1  -  5). 
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SUMMARY 


Two  novel  strains  of  Propionigenium  maris  able  to  reductively  debrominate  2,4,6- 
tribromophenol  (TBP)  to  monobromophenols  were  isolated  from  marine  hemichordate  and 
polychaete  burrows.  These  two  strains,  DSL-1  and  ML-1,  were  anaerobic,  nonmotile  rods  that 
stained  Gram  negative  and  required  0.05%  yeast  extract  for  growth.  Strain  DSL-1  fermented 
pyruvate  and  succinate  to  butyrate  and  strain  ML-1  fermented  glucose  and  succinate  primarily  to 
propionate.  No  inorganic  terminal  electron  acceptors  were  identified.  The  pH  and  temperature 
optima  for  growth  of  strain  DSL-1  were  7.6  and  30°  C;  7.0  and  32°  C  for  strain  ML-1  with 
doubling  times  of  0.32  h  and  0.30  h,  respectively.  Both  strains  required  2-3%  (w/v)  NaCl  for 
optimal  growth.  Morphological  and  physiological  features,  as  well  as  the  results  of  16S  rDNA 
sequence  analysis,  showed  these  to  be  new  strains  of  Propionigenium  maris.  Because  they  differ 
from  the  P.  maris  type  strain  (DSM  9537)  in  a  number  of  respects,  including  their  ability  to 
rapidly  debrominate  di-  and  tribromophenols,  and  in  their  specific  habitats,  the  species 
description  is  amended  to  include  these  ecologically  important  properties. 
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INTRODUCTION 


Biogenic  bromoaromatic  compounds  are  broadly  distributed  in  intertidal  marine 
sediments  (King,  1986;  D.  E.  Lincoln,  K.  T.  Fielman,  R.  L.  Marinelli,  and  S.  A.  Woodin, 
unpublished  data;  Steward  et  al,  1992;  Woodin  et  al,  1987).  A  variety  of  infaunal  hemichordate 
and  polychaete  worms  produce  bromoaromatic  compounds,  including  but  not  limited  to 
bromophenols,  bromopyrroles,  bromoindoles,  and  bromobenzylalcohols  (Fielman  et  ah,  1999; 
Woodin  et  al. ,  1987).  These  bromometabolites  are  produced  through  the  action  of 
haloperoxidases  (Chen  et  al.  ,1991)  and  reach  very  high  levels  in  the  worm  tissues  (Fielman  & 
Target!,  1995;  Yoon  et  al.,  1994).  The  highest  concentrations  of  bromometabolites  external  to 
the  animals  are  found  in  sediments  lining  the  animals'  burrows  (King,  1986;  D.  E.  Lincoln,  K.  T. 
Fielman,  R.  L.  Marinelli,  and  S.  A.  Woodin,  unpublished  data).  Turnover  of  bromophenols  in 
shallow  marine  sediments  can  be  rapid  (Steward  &  Lovell,  1997)  and  the  burrow  linings,  due  to 
the  abundance  of  carbon  cosubstrates  and  potential  electron  acceptors,  and  high  levels  of 
bacterial  biomass  (Steward  et  al.,  1996;  Phillips  &  Lovell,  1999),  are  expected  to  be  sites  of 
enhanced  degradation  of  bromoaromatic  compounds.  The  microbiota  inhabiting  these  burrow 
lining  biofilms  are  thus  of  considerable  interest  for  detailed  studies  of  degradation  processes  for 
haloaromatic  compounds.  Presently  little  is  known  about  the  diversity  of  species  carrying  out 
reductive  dehalogenation  processes  in  marine  systems. 


The  large  extent  of  contaminated  sediments  and  the  diversity  of  microenvironments 
available  in  coastal  marine  systems  implies  the  existence  of  a  plethora  of  haloaromatic  degrading 
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species  there.  However,  although  mixed  culture  and  sediment  slurry  studies  of  marine 
haloaromatic  degradation  are  abundantly  documented  in  the  literature  (e.g.  Abrahamsson  & 
Klick,  1991;  Allard  eta/.,  1991;  Boothe  et  a/.,  1997;  Genthner  e/a/.,  1989;  Haggblom &Young, 
1990,  1995;  King,  1988;  Monserrate  &  Haggblom,  1997;  Remberger  et  al,  1986),  and  many 
important  degradative  reactions  are  known  (see  El  Fantroussi  et  al,  1998;  Fetzner  &  Lingens, 
1994;  Haggblom,  1992;  Hale  et  al,  1994;  Janssen  et  al,  1994;  Mohn  &  Tiedje,  1992  for  recent 
reviews),  marine  strains  highly  active  against  multiply  substituted  haloaromatics  have  proyen 
difficult  to  enrich  and  isolate  into  pure  culture.  This  is  particularly  true  of  reductively 
dehalogenating  organisms,  which  can  remove  halogen  atoms  from  the  aromatic  ring;  This 
process  is  crucial  in  the  consortial  degradation  of  multiply  substituted  haloaromatics  since 
compounds  with  few  substituent  halogen  atoms,  particularly  monohalogenated  aromatics,  are 
mineralized  much  more  readily  than  multiply  substituted  compounds  with  the  same  aromatic  ring 
structure  (Neilson,  1990,  Utkin  et  al,  1995).  No  pure  culture  able  to  dechlorinate  multiply 
chlorinated  chloroaromatic  compounds  has  been  reported  from  marine  systems  and  only  one 
marine  strain  able  to  debrominate  multiply  brominated  aromatics  has  been  described  to  date 
(Steward  et  al,  1995).  This  limited  roster  of  dehalogenators  restricts  studies  of  the  physiology 
and  biochemistry  of  reductive  dehalogenation  in  marine  microorganisms,  a  reaction  catalyzed  by 
several  classes  of  enzymes  differing  in  structure  and  in  some  cases  in  reaction  mechanism  (Y.  P. 
Chen,  J.  Watson ,  and  C.  R.  Lovell,  unpublished  data;  Holliger  et  al,  1999,  Ldffler  et  al.,  1996; 
Ni  et  al,  1995;  Steward  et  al,  1995).  Here  we  present  the  characterization  and  phylogenetic 
affiliation  of  two  reductively  debrominating  marine  anaerobe  strains  of  Propionigenium  mans. 
We  amend  the  description  of  this  species  since  the  two  novel  strains  are  capable  of  rapid 
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reductive  debromination  of  tri-  and  dibromophenols,  producing  monobromophenols  as  terminal 
products,  a  property  not  described  for  the  type  strain. 

MATERIALS  AND  METHODS 

Materials  and  strains.  Bromophenol  substrates  and  standards  were  obtained  from  Aldrich 
Chemical  and  used  without  further  purification.  Propionigenium  mans  (DSM  9537)  was 
obtained  from  the  Deutsche  Sammlung  von  Mikroorganismen.  The  isolation  of  strain  DSL-1  has 
been  described  elsewhere  (Steward  et  al,  1995).  Samples  for  enrichment  of  strain  ML-1  were 
collected  from  the  North  Inlet  salt  marsh,  Georgetown,  South  Carolina,  USA.  (30°20'N, 

79°1 0'W).  These  samples  consisted  of  burrow  linings  of  the  capitellid  polychaete  Notomastus 
lobatus,  a  common  inhabitant  of  clay  rich  intertidal  mud  flats  in  this  system  and  a  well-studied 
producer  of  2,4,6-tribromophenol,  2,4-dibromophenol,  and  4-bromophenol  (Chen  et  a/.,  1991; 
Lovell  et  al,  1999;  Steward  et  al,  1992, 1996;  Steward  &  Lovell,  1997;  Yoon  et  al,  1994).  The 
burrow  scrapings  were  transported  on  ice  to  Columbia,  South  Carolina,  where  enrichment 
cultures  were  prepared  following  previously  described  methods  (Steward  et  al,  1995). 
Enrichment  media  had  pH  values  ranging  from  5.5  to  7.5,  and  all  contained  1  mM  2,4,6- 
tribromophenol  (TBP);  a  concentration  comparable  to  that  of  total  bromophenols  in  N.  lobatus 
burrow  linings  (D.  E.  Lincoln,  K.  T.  Fielman,  R.  L.  Marinelli,  and  S.  A.  Woodin,  unpublished 
data).  These  enrichment  media  also  contained  acetate,  lactate,  sodium  benzoate,  or  succinate  (40 
mM  final  concentration)  as  supplemental  carbon  source.  The  cultures  were  incubated  in  the  dark 
at  32°  C.  Gro'wth  was  seen  within  a  few  days  in  all  cultures  except  those  containing  sodium 
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benzoate.  All  growing  cultures  actively  debrominated  TBP  as  assessed  using  gas 
chromatography  (Stewart  et  al,  1995).  The  pH  7.0,  succinate  supplemented  culture  had  few 
microscopically  distinguishable  cell  types  and  was  selected  for  pure  culture  isolation  of 
debrominators  using  the  agar  shake  method  (Steward  et  al. ,  1995).  Colonies  were  visible  within 
one  week  and  were  aseptically  transferred  into  fresh  succinate  medium.  This  process  was 
repeated  four  times,  with  each  successive  culture  checked  for  homogeneity  by  phase  contrast 
microscopy.  All  cultures  were  homogeneous  with  respect  to  cell  morphology  and  were 
apparently  identical  to  each  other.  One  pure  culture  from  this  group  was  selected  arbitrarily  and 
designated  strain  ML- 1 . 

Growth  media.  DSL-1 :  The  growth  medium  was  a  minor  modification  of  the  recipe  of  Steward 
et  al.  (1995)  and  contained  the  following:  (g  •  1’)  KH2PO4, 0.2;  NH4CI,  0.3;  KCl,  0.5;  CaCl2  • 
2H2O,  0.15;  NaCl,  20;  MgClj  •  6H2O,  3.0;  Bacto  yeast  extract,  0.5;  SL-9  trace  metal  solution 
(Tschech  and  Pfennig,  1984),  1  ml  •  1''.  The  basal  medium  was  adjusted  to  pH  5.5  and,  after 
autoclaving,  the  medium  was  supplemented  with  sterile  filtered  1  mM  2,4,6-tribromophenol, 
autoclaved  cysteine  sulfide  (0.0028%  cysteine  -  0.0028%  sodium  sulfide,  w/v,  final 
concentration),  and  autoclaved  NaHCOj  (50  mM  final  concentration).  In  some  experiments,  a 
supplemental  carbon  source  (acetate,  benzoate,  citrate,  crotonate,  ethanol,  fiimarate,  glucose, 
lactate,  malate,  oxalate,  propionate,  pyruvate,  or  succinate  at  4  -  20  mM  final  concentration) 
and/or  potential  terminal  electron  acceptor  (nitrate,  fiimarate,  and  sodium  thiosulfate  at  5  mM 
final  concentration)  were  added.  Routine  cultivation  of  strain  DSL-1  was  performed  using  20 
mM  succinate  as  supplemental  carbon  source  with  no  added  terminal  electron  acceptor.  The 
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complete  medium  was  adjusted  to  pH  7.5. 

ML-1:  The  growth  medium  for  strain  ML-1  was  derived  from  the  DSL-1  medium  recipe, 
substituting  4-(2-hydroxyethyl)-l-piperazineethanesulfonic  acid  (HEPES,  4.8  g  •  T’)  for  the 
carbonate  buffer,  and  adding  Bacto  Casitone  (1.0  g  •  T')  as  an  additional  supplement.  Routine 
cultivation  employed  the  ML-1  basal  medium  supplemented  with  40  mM  succinate.  Culture 
growth  was  monitored  by  direct  microscopic  counting  using  a  Petroff-Hauser  counting  chamber. 

pH  and  temperature  ranges  and  NaCI  requirements.  The  pH  and  temperature  ranges  for 
growth  of  strains  DSL-1  and  ML-1  were  determined  in  triplicate  for  cultures  growing  on 
succinate  supplemented  basal  media.  Optimum  pH  and  NaCl  requirement  determinations  were 
performed  at  room  temperature  for  strain  DSL-1  and  at  32°  C  for  strain  ML-1 .  The  pH  range 
tested  for  strain  DSL-1  was  from  6.5  through  8.5  and  for  strain  ML-1  from  5.5  through  7.5.  The 
temperature  range  for  growth  of  strain  DSL-1  was  determined  at  pH  7.5  over  a  range  from  10°  C 
through  50°  C.  The  temperature  range  for  strain  ML-1  was  determined  at  pH  7.0  and 
temperatures  between  20°  C  and  40°  C  were  used.  Both  strains  were  tested  for  growth  within  the 
ranges  specified  at  five  degree  temperature  intervals.  Once  the  approximate  temperature  optima 
were  determined,  the  strains  were  retested  using  two  degree  temperature  intervals  near  the 
optima.  The  NaGl  requirement  for  each  organism  was  determined  by  varying  the  NaCl  content 
of  the  succinate  supplemented  basal  medium  recipes.  The  unmodified  media  both  contained  2% 
NaCl  (w/v).  We  inoculated  from  these  2%  NaCl  media  into  media  containing  different  NaCl 
levels  from  0  to  6%. 
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Determination  of  fermentation  products.  Fermentation  products  were  determined  using  a 
Hewlett  Packard  1100  series  high  performance  liquid  chromatography  system  equipped  with  a  25 
cm  Microsorb-MV  Cl 8  reversed-phase  column  (RaininA^arian)  and  using  a  0.01  N  H2SO4 
mobile  phase  at  a  flow  rate  of  1  ml  •  min'.  Detection  was  by  absorbance  at  210  nm.  Retention 
times  of  products  were  compared  to  those  of  known  standard  compounds,  including  acetate, 
propionate,  and  butyrate.  Compound  quantification  was  based  on  absorbances  of  known 
quantities  of  authentic  compounds.  Fermentation  products  were  determined  during  early  and  late 
exponential  phase,  again  during  early  stationary  phase,  and  finally  late  in  stationary  phase. 

Dehalogenation  assays.  Dehalogenation  was  determined  at  timed  intervals  throughout  growth 
by  measuring  dehalogenation  products  by  HPLC.  Samples  (1  ml)  were  taken  and  mixed  with  an 
equal  volume  of  100%  acetonitrile.  After  centrifugation  to  remove  precipitated  proteins  and 
polysaccharides,  the  samples  were  assayed  for  bromophenols  using  the  Hewlett  Packard  1 1 00 
series  HPLC  system  with  a  1 5  cm  HP  Cl 8  reversed-phase  column  and  using  a 
methanoLHjOiacetic  acid  (50:49:1)  mobile  phase  at  a  flow  rate  of  1.5  ml  •  min  '.  Detection  was 
measured  by  absorbance  at  280  nm.  Retention  times  of  dehalogenation  products  were  compared 
to  those  of  authentic  standard  compounds,  including  2-bromophenol,  4-bromophenol,  2,4- 
dibromophenol,  2,6-dibromophenol,  and  2,4,6-tribromophenol  (TBP).  Compound  quantification 
was  based  on  absorbances  of  known  quantities  of  authentic  compounds. 

Determination  of  mol%  G  +  C  content.  DNA  was  isolated  using  a  modification  of  the  Marmur 
(1961)  procedure,  adding  a  cetyltrimethylammmonium  bromide  extraction  step  for  removal  of 
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contaminating  polysaccharides  and  proteins  (Ausubel  et  al,  1987).  The  mol%  G  +  C  contents  of 
each  strain  were  determined  after  PI  nuclease  digestion  and  HPLC  separation  of  nucleosides 
(Mesbah  et  al. ,  1989;  Whitman  et  al. ,  1 986). 

16S  rDNA  sequence  determination  and  analysis.  Genomic  DNA  was  extracted  and  the  gene 
coding  for  16S  rRNA  (16S  rDNA)  was  amplified  as  described  previously  (Rainey  et  al,  1996). 
Purified  PCR  products  were  directly  sequenced  using  the  Taq  DyeDeoxy  Terminator  Cycle 
Sequencing  Kit  (Applied  Biosystems).  Sequence  reactions  were  electrophoresed  using  Applied 
Biosystems  models  373 A  and  310  DNA  sequencers.  BLAST  analysis  was  preformed  on  the 
sequence  data  to  determine  the  phylogenetic  grouping  to  which  the  strains  were  most  closely 
related  (Altschul  et  al.,  1990).  The  16S  rDNA  sequences  of  strains  DSL-1  and  ML-1  were 
manually  aligned  with  sequences  of  representatives  of  the  low  G+C  Gram  positive  subphylum 
using  the  ae2  editor  and  sequence  similarity  values  calculated  by  pairwise  sequence  comparison 
(Maidak  et  al.,  1994).  The  16S  rDNA  sequences  determined  in  this  study  are  available  from 
EMBL  under  the  accession  numbers  Y16799  (strain  DSL-1)  and  Y16800  (strain  ML-1). 

RESULTS  AND  DISCUSSION 


Phylogenetic  analyses.  16S  rDNA  sequences  comprising  1441  nucleotides  between  positions  31 
and  1493  {Escherichia  coli  numbering)  were  determined  for  strains  DSL-1  and  ML-1 . 
Comparative  analysis  of  the  16S  rDNA  sequences  with  those  sequences  available  for 
representatives  of  the  major  bacterial  groups  using  BLAST  analyses  indicated  that  both  strains 
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were  members  of  the  low  G+C  Gram  type  positive  subphylum.  Alignment  of  the  two  sequences 
with  those  available  for  the  members  of  the  low  G+C  Gram  type  positive  subphylum  and 
analysis  of  all  nucleotide  positions  showed  them  to  be  most  closely  related  to  members  of  the 
genus  Propionigenium.  The  16S  rDNA  sequence  of  strain  DSL-1  was  found  to  have  99.7% 
similarity  to  that  of  strain  ML-1 .  The  16S  rDNA  sequence  similarity  between  strain  DSL-1  and 
P.  maris  was  99.6,%  while  that  to  P.  modestum  was  96.0%.  The  16S  rDNA  sequence  of  strain 
ML-1  had  99.9%  similarity  to  P.  maris  and  95.7%  to  P.  modestum. 

Characterization  of  strains  DSL-1  and  ML-1. 

Cell  morphology.  Cells  of  strain  DSL-1  in  exponential  growth  phase  were  short  rods  averaging 
1  pm  in  length  and  0.5  pm  in  diameter  (Table  1).  Stationary  phase  cells  were  shorter,  some 
almost  spherical  (Figure  1).  This  strain  displayed  asymmetric  cell  division,  with  one  daughter 
cell  noticeably  larger  than  the  other  during  stationary  phase  (Figure  1).  Strain  ML-1  was  also  rod 
shaped  and  0.5  pm  in  diameter,  but  cell  length  during  exponential  growth  was  variable,  ranging 
from  2  to  10  pm.  Both  strains  became  shorter  and  more  uniform  during  stationary  phase,  with 
ML-1  length  averaging  2  pm  in  length.  ML-1  did  not  display  asymmetric  cell  division. 
Sporulation  was  not  observed  in  either  strain  and  is  not  a  property  of  the  genus  Propionigenium 
(Schink  &  Pfennig,  1982).  Strain  ML-1  contained  dense  inclusions  that  were  most  numerous  in 
stationary  phase  cells  (Figure  2).  No  such  inclusions  were  found  in  DSL-1  in  any  growth  phase. 

Gram  stain  reaction.  The  Gram  stain  reaction  was  negative  but  electron  microscopy  of  osmium 
tetraoxide  stained  thin  sections  revealed  a  single  dark-staining  cell  envelope  layer  (data  not 
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shown),  consistent  with  a  Gram  type  positive  cell  wall  structure.  No  lipopolysaccharide  (LPS), 
the  characteristic  compound  of  Gram  type  negative  cell  walls,  was  detected  using  the  LPS- 
polymyxin  B  assay  of  Wiegel  &  Quandt  (1982).  Other  Propionigenium  strains  also  stain  Gram 
negative  (Schink  &  Pfennig,  1982). 

Growth  properties.  Both  strains  were  obligately  anaerobic.  The  pH  range  for  growth  of  strain 
DSL-1  was  6.5  to  8.5,  with  optimal  growth  at  pH  7. 5-7. 8.  The  pH  range  supporting  growth  of 
strain  ML-1  was  5.5  to  7.5,  and  the  optimum  pH  was  7.0.  The  lower  pH  optimum  for  strain  ML- 
1  may  reflect  the  stagnant  conditions  in  N.  lobatus  burrows  (Steward  et  al,  1996).  Burrow 
stagnation  permits  accumulation  of  volatile  fatty  acids  (chiefly  propionate,  lactate,  and  acetate) 
produced  by  the  worm  under  conditions  of  anaerobiosis  (Schdttler  et  a/.,  1983;  Surholt,  1977) 
and  consequently,  substantial  decreases  in  burrow  water  pH  (Scott  1976).  The  temperature  range 
for  growth  of  strain  DSL-1  was  10  to  35°  C,  and  the  optimum  temperature  was  30°  C.  The 
optimum  temperature  for  growth  of  strain  ML-1  was  similar,  32°  C.  The  shortest  doubling  times 
observed  for  strains  DSL-1  and  ML-1  at  optimum  pH  and  growth  temperature  were  0.30  -  0.32  h. 
Strain  DSL-1  grew  well  within  a  NaCl  concentration  range  of  1-5%  with  optimal  grov^h  at  2- 
3%  NaCl  and  no  growth  at  0  or  6%.  Strain  ML-1  grew  well  from  1-4%  NaCl  with  optimal 
growth  at  3%  and  no  growth  at  0  or  5%.  These  results  are  consistent  with  the  salinities  observed 
in  the  euhaline  saltmarsh  environment  from  which  these  strains  were  isolated.  The  only  defined 
carbon  sources  that  were  utilized  in  the  presence  of  0.05%  yeast  extract  were  pyruvate  and 
succinate  for  strain  DSL-1  and  glucose  and  succinate  for  strain  ML-1.  The  growth  rate  of  strain 
DSL-1  was  increased  over  that  in  the  basal  medium  by  addition  of  succinate  or  pyruvate  at 
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concentrations  up  to  20  mM,  but  not  above,  while  the  growth  rate  of  strain  ML-1  was  increased 
by  succinate  at  concentrations  up  to  40  mM.  Substrates  that  did  not  measurably  increase  the 
growth  rate  of  either  strain  in  the  presence  of  0.05%  yeast  extract  included  acetate,  benzoate, 
citrate,  crotonate,  ethanol,  fumarate,  lactate,  malate,  oxalate,  and  propionate.  In  addition,  strain 
DSL-1  did  not  utilize  glucose  and  ML-1  did  not  utilize  pyruvate.  Strain  DSL-1  growing  on 
either  pyruvate  or  succinate  showed  exponential  decline  in  cell  numbers  in  stationary  phase.  In 
the  absence  of  these  supplemental  carbon  sources  the  decline  in  stationary  phase  cell  numbers 
was  less  marked.  A  similar  effect  was  observed  for  glucose-grown  strain  ML-1,  which  goes  into 
exponential  decline  after  stationary  phase.  These  stationary  phase  cultures  could  not  always  be 
successfully  subcultured.  Strain  ML-1  cultures  grown  with  succinate  could  be  successfully 
subcultured  after  as  long  as  4  months  if  they  were  kept  at  room  temperature.  No  appreciable 
change  in  growth  rate  or  yield  of  either  strain  resulted  from  addition  of  nitrate,  fumarate,  or 
thiosulfate  as  potential  terminal  electron  acceptors.  The  fermentation  product  of  strain  DSL-1  at 
all  growth  stages  was  butyrate.  Strain  ML-1  consistently  produced  propionate.  Similar  to  strain 
ML-1,  the  P.  marts  type  strain  (DSM  9537)  utilizes  succinate  and  glucose  (as  well  as  several 
other  substrates),  requires  yeast  extract,  and  consistently  produces  propionate  as  its  fermentation 
product  (Janssen  &  Liesack,  1995). 

Debromination.  Tribromophenol  was  rapidly  debrominated  by  both  strains  (Steward  et  al, 

1995  and  Figure  3).  DSL-1  produced  primarily  2-bromophenol  as  the  terminal  debromination 
product;  ML-1  produced  a  mixture  of  2-bromophenol  and  4-bromophenol.  The  discrepancy 
between  the  previous  report  of  4-bromophenol  production  by  strain  DSL-1  (Steward  et  al,  1995) 
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and  the  current  observation  of  2-bromophenol  is  difficult  to  explain  since  in  both  cases  the 
product  identities  were  confirmed  by  mass  spectroscopy.  This  change  may  be  due  to  a  change  in 
the  culture  between  the  initial  determinations  and  the  present  work.  Neither  strain  could  utilize 
2-  or  4-bromophenol  as  carbon  substrates.  3-bromophenol  was  not  tested.  Neither  strain 
dechlorinated  chlorophenols.  Similar  substrate  spectra,  i.e.  the  capacity  to  debrominate  but  not 
to  dechlorinate  the  corresponding  chlorinated  compounds,  has  been  observed  with  2,4,6- 
tribromobiphenyl  dehalogenating  enrichment  cultures  (Chuang,  1997).  However,  the 
corresponding  2,4,6-trichlorobiphenyl  dehalogenating  enrichment  could  dehalogenated  both  the 
chlorinated  and  structurally  analogous  brominated  compounds  (Chuang,  1997),  indicating  the 
presence  of  at  least  two  different  assemblages  of  dehalogenating  bacteria,  of  which  one  could 
only  debrominate. 

Bromophenols  supplied  at  millimolar  concentrations  did  not  affect  yield  or  doubling  time.  The 
absence  of  growth  stimulation  or  change  in  growth  yields  indicates  that  debromination  by  these 
strains  does  not  lead  to  energy  production.  This  is  in  contrast  to  the  dehalogenation  of 
chloro(bromo)phenols  by  the  Gram  type  positive,  sulfite  reducing  Desulfitobacterium 
dehalogenans,  which  was  isolated  from  freshwater  sediments,  and  in  which  dehalogenation  leads 
to  1  ATP  formed  for  each  reductively  removed  halogen  atom  and  thus  to  increased  growth  yields 
in  comparison  to  the  fermentation  of  pyruvate.  The  energy  is  presumably  coupled  to  electron 
transport  phosphorylation.  This  mechanism  enables  this  dehalogenating  bacterium  to  grow  on 
hydrogen  and  formate  while  using  the  halogenated  phenols  as  electron  acceptors  (Mackiewicz  & 
Wiegel,  1998).  The  debromination  of  2,4-dibromobiphenyl  also  supports  enhanced  growth  of 
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reductively  dehalogenating  microorganisms  in  sediment  samples  (Wu  et  al,  1999).  P.  marts, 
the  closest  related  organism  to  the  two  novel  isolates  based  on  16S  rDNA  sequence  analysis,  was 
also  tested  for  dehalogenation  activity.  Growth  of  unadapted  P.  marts  type  strain  DSM  9537  was 
strongly  inhibited  by  1  mM  TBP.  After  adaptation  through  successive  subculturing  using  media 
containing  1 00  pM  through  500  pM  TBP,  type  strain  DSM  9537  could  be  grown  in  the  presence 
of  up  to  500  pM  TBP  and  trace  amounts  of  the  debromination  terminal  product,  2-bromophenol, 
were  detected.  No  indication  of  growth  stimulation  or  evidence  that  debromination  is  coupled  to 
energy  production  were  obtained  for  any  P.  marts  strains.  Proptontgentum  species  are  regarded 
as  purely  fermentative  anaerobes  (Schink  &  Pfennig,  1982;  Janssen  &  Liesack,  1995).  No 
increases  in  cell  yields  by  P.  modestum  (Schink  &  Pfennig  1982)  or  any  strain  of  P.  marts  in 
response  to  added  nitrate,  fumarate,  or  thiosulfate  has  been  observed.  This  indicates  that  no 
energy  yielding  electron  transport  phosphorylation  system  is  functioning  in  these  species. 

G+C  content.  The  G+C  content  of  genomic  DNA  from  DSL-1  is  39  ±  0.3  mol%  (SD,  n  =  6),  41 
±  0.1  mol%  (SD,  n  =  4)  for  ML-1  DNA.  These  values  are  in  good  agreement  with  G+C  contents 
of  other  Proptontgentum  strains  (Table  1). 

In  conclusion,  the  high  16S  rDNA  sequence  similarity  found  between  strains  DSL-1  and  ML-1 
and  the  two  described  species  of  the  genus  Proptontgentum  clearly  indicate  their  membership  in 
the  genus  Proptontgentum.  Based  on  similarities  to  the  P.  marts  type  strain  and  phylogenetic 
analysis  we  currently  place  the  new  isolates  as  novel  strains  in  the  species  Proptontgentum 
marts.  The  new  strains  have  a  high  degree  of  16S  rDNA  sequence  similarity  to  that  of  the  P. 
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marts  type  strain  (99.6-99.9%),  but  differ  from  the  latter  organism  in  several  physiological 
features  including  cell  size,  substrate  preferences,  temperature  and  pH  optima  and,  most 
importantly,  ability  to  rapidly  debrominate  TBP.  The  infaunal  burrow,  a  dynamic 
microenvironment  that  supports  high  bacterial  biomass  and  activities,  may  promote  rapid 
evolution  of  these  Propionigenium  strains,  resulting  in  substantial  physiological  variability 
without  concurrent  changes  in  16S  rDNA  sequence.  Our  identification  of  these  isolates  as  strains 
of  P.  marts  is  conservative  and  based  in  large  measure  on  16S  rDNA  sequence  homology. 

The  new  P.  marts  strains  have  been  isolated  from  infaunal  burrows  with  important  differences  in 
several  microenvironmental  features.  The  size,  shape,  and  placement  in  the  sediments  of  the 
burrows,  the  texture  and  organic  matter  content  of  sediment  in  the  burrow  linings,  and  the 
porosity  of  the  surrounding  sediments  (which  affects  solute  transport  rates)  all  differ  (see 
Steward  et  al,  1996  for  further  discussion).  The  feeding  modes  of  the  host  invertebrates,  their 
irrigation  behaviors,  and,  likely,  the  availability  of  organic  substrates  released  by  the  hosts  differ 
as  well.  Given  the  great  variety  of  burrowing  marine  invertebrate  taxa,  the  many  types  of  burrow 
structures  they  build,  and  their  range  of  microbiologically  significant  behaviors  (burrow 
maintenance  activities,  irrigation  regimes,  feeding  strategies,  etc.)  it  is  probable  that  many 
distinct  strains  of  P.  marts,  and  likely  new  Proptontgentum  species,  occur  in  marine  infaunal 
burrows.  Burrows  of  bromometabolite  producing  infauna  may  also  harbor  additional 
dehalogenating  strains  and  species.  Further  work  is  needed  for  evaluation  of  the  importance  of 
this  group  of  organisms. 
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Emendation.  The  species  description  for  P.  maris  (Janssen  &  Liesack,  1995)  is  emended  for  the 
property  of  reductive  dehalogenation  of  bromophenols  found  in  the  linings  of  some  marine 
infaunal  burrows.  In  contrast  to  strains  isolated  from  the  burrows,  which  are  regarded  as  a 
special  habitat  for  some  P.  maris  strains,  the  type  strain  DSM  9537  exhibits  only  minor 
debromination  activity  and  only  after  prolonged  adaptation.  Furthermore,  some  strains  may  not 
utilize  sugars  and  butyrate  can  be  the  major  fermentation  product  from  pyruvate  and  succinate, 
although  the  type  strain  forms  propionate  as  the  major  product  from  glucose  and  succinate. 
Depending  on  the  isolate,  the  cell  size  varies  between  0.5  and  1  pm  in  diameter  and  1  and  10  pm 
in  length. 
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Table  1.  Selected  properties  of  Propionigenium  marls  strains. 
P.  modestum  is  included  for  comparison. 


Property 

P.  marls 
DSL-1 

P.  marls 
ML-1 

P.  marls* 
DSM  9537 

P.  modestum^ 

Cell  size  (Hm)* 

1  X  0.5 

2  X  0.5  ■ 

2 . 5  X  1 

2  X  0 . 5 

Inclusions 

none 

present 

none 

.  none 

Carbon  substrates 

succinate 

pyruvate 

succinate 

glucose 

succinate 

pyruvate 

fumarate 

citrate 

3-hydroxy- 

butyrate 

glucose 

fructose 

maltose 

succinate : 
pyruvate 
fumarate 
malate 

oxaloacetate 

L-aspartate 

Fermentation 

Products^ 

butyrate 

propionate 

propionate 

(butyrate) 

(acetate) 

(lactate) 

(formate) 

(ethanol) 

propionate 

(acetate) 

Optimum  pH 

7. 5-7. 8 

7.0 

6. 9-7. 7 

7. 1-7.7 

Optimtum 

temperature  {  C) 

:  28 

32 

34-37 

33 

Mol%  G+C 

39 

41 

4  0 

:.33.9 

Data  from  reference  Janssen  &  Liesack  (1995) . 

-Data  from  reference  Schink  &  Pfennig  (1982)  . 

*  Cell  dimensions  (length  x  width)  for  exponentially  growing 
cells.  All  Propionigenium  cells  are  rods 

§  Production  of  butyrate,  acetate,  lactate,  formate,  and  ethanol 
by  P.  marls  and  acetate  by  P.  modestum  is  dependent  bn  the  carbon 
source . 
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FIGURE  LEGENDS 


Figure  1.  Phase  contrast  photomicrograph  of  Propionigenium  maris  strain  DSL-1  during 
stationary  phase.  Note  asymmetric  cell  division  (arrows). 

Figure  2.  Phase  contrast  photomicrograph  of  Propionigenium  maris  strain  ML-1  during 
stationary  phase.  Note  inclusions  (arrows). 

Figure  3.  Growth  (A)  and  debromination  of  2,4,6-tribromophenol  (B)  by  Propionigenium  maris 
strain  ML-1 .  (Open  circles)  2,4,6-tribromophenol,  (closed  squares)  2,4-dibromophenol,  (open 
triangles)  4-bromophenol,  (closed  triangles)  2-bromophenol.  Values  are  means  of  three  cultures 
incubated  in  parallel.  Error  bars  represent  one  standard  deviation.  For  symbols  with  no  error 
bars,  error  bars  are  smaller  than  the  symbols. 
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